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THE LIMITING SENSITIVITY OF 
SEISMIC DETECTORS* 


ALFRED WOLFt 


ABSTRACT 


The theory of Brownian motion is applied to the problem of determination of the 
limiting sensitivity of seismic detectors. It is shown that a certain minimum suspended 
mass is requird for the recording of small ground motions. Electromagnetic geophones 
are studied in detail, and it is shown how spontaneous thermal fluctuations in voltage 
limit the performance Of these instruments. Under certain conditions, a geophone may 
be treated as a power generator, and the necessary suspended mass is then determined 
by the power requirements of the recording apparatus. 


The sensitivity of instruments for the measurement of very small 
physical quantities is often limited by irregular thermal motion, the 
Brownian motion, of some of their component parts. The case of gal- 
vanometers is the best known, and has been studied in detail both 
experimentally and theoretically.! Similarly, the minimum amplitude 
of ground motion which can be detected by a mechanical seismograph 
or an electromagnetic geophone of a given design is limited by the 
Brownian motion of the mass which is suspended in the instrument. 

The purpose of this article is the application of the theory of 
Brownian motion and of other similar phenomena to the principles of 
design of seismic detectors, more particularly those employed in geo- 
physical exploration work. It will be shown that the limiting sensi- 
tivity of most geophones employed for the recording of reflections is 
so great that no difficulty should be experienced in detecting motions 
as small as ground unrest even in quiet locations. Such geophones are, 
however, designed to respond only to frequencies of the order of 30 
cycles/sec. or higher. The problem of detecting small ground motions 
becomes much more difficult when the frequency is of the order of 1 
to 5 cycles/sec., as it may be in long distance refraction work. It is 


* Presented at the Annual Meeting, Houston, Texas, April 2, 1941. 
} Geophysical Research Corporation, Tulsa, Oklahoma. 
1 Cf. Barnes and Silverman: Reviews of Modern Physics, July, 1934. 
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in this region that the theory to be presented would be of some value. 

The main results, as applied to the usual type of geophone, may 
be stated as follows: The power which a geophone derives from ground 
motion is proportional to the product of the suspended mass of the 
geophone and the square of the amplitude of ground motion; it is, 
incidentally, also a function of frequency. The average energy of ther- 
mal motion of the geophone mass is a function only of the temperature 
and of the frequency range. It follows that the mass of a geophone 
must exceed a certain minimum value if a given ground motion is to 
be registered. 

In the first section of this article, the limiting sensitivity of seismic 
detectors of the displacement and of the velocity types will be deduced 
from very general principles. In the second section, the particular case 
of an electromagnetic geophone will be considered from a more prac- 
tical point of view. It will be shown that thermal fluctuation voltage 
is one of the important factors which limit the sensitivity of this type 
of geophone. In the third and concluding section, the geophone is 
treated as a generator of power which must be furnished the recording 
apparatus for satisfactory performance. 

Only instruments which are critically damped will be discussed; 
these are the most commonly employed for the registration of earth- 
quakes and in geophysical exploration work. 


THE BROWNIAN MOTION OF THE SUSPENDED 
MASS OF A SEISMIC DETECTOR 


The equation of motion of the suspended mass of a critically 
damped seismic detector can be written 


+ + wo2x = F/M (1) 


where x is the displacement relative to the ground, M is the mass, 
wo/2m=fo is the natural frequency of the seismic detector. F is the 
force due to ground motion and to other causes; in computations of 
_ Brownian motion the force F may, for instance, be the result of mo- 
lecular motion. 

When the force F is due to harmonic ground motion of amplitude 
A and frequency w/2r=f 


F = (2) 
and it follows from (1) 
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and further 


I wtA? (4) 
x = 
2 (w? + w?)? 
and 
I 
(5) 


(w? + 


x? and #? are respectively the mean square values of the displacement 
and of the velocity of the suspended mass of the seismic detector. 

The irregular thermal motion of the suspended system of the seis- 
mic detector can be dealt with according to well-known principles of 
statistical mechanics. In this case, when the motion is governed by (1), 
the statement can be made that the mean kinetic energy and the mean 
potential energy of the suspended mass, which are due to thermal 
motions, are functions of the absolute temperature of the instrument 
and can be written as 


= = RT (6) 


Xp, Xp are respectively the Brownian motion displacement and ve- 
locity; k is a universal constant (k=1.38X10~"* erg/deg K), T is the 
absolute temperature in degrees Kelvin. 

If the motion of the suspended mass of the seismic detector were 
recorded without frequency discrimination, it would now be a simple 
matter to write out the relation for the limiting sensitivity of the 
seismograph. For, it would only be necessary to state that the mean 
square motion of the suspended mass, as given by (4) or (5), must 
exceed the mean square Brownian motion by a factor a which is de- 
rived from experience. Practically, however, all types of recorders 
have some frequency limitation, so that the recorded Brownian motion 
appears to be smaller than that given by (6). It will be assumed in the 
following that the frequency discrimination of the recording device 
is equivalent to that of a low-pass filter with a cut-off frequency 
2a = fi. 

The analytic expression for the statement that the mean square 
motion due to ground acceleration must exceed the mean square 
Brownian motion by at least the factor a can be written 


2kT 2 2 2 
(- ) (7) 


2 w2 


Wo 
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for a displacement type detector, or 


kT 2 2 
(=) a (7’) 


w? w? 


for velocity recording. 

The factors G; and G2 are introduced to take account of the fre- 
quency cut-off of the recorder; they are functions of w; and wo of the 
order of unity, as will be demonstrated. 

Experience with galvanometers indicates that a motion of ampli- 
tude four or five times greater than the amplitude of Brownian motion 
can be registered with certainty; it will, accordingly, be assumed in 
numerical estimates of sensitivity that 

a= 20 
since it is proportional to the square of the amplitude. 

The equality sign in (7) and (7’) is used when A is identified with 
the minimum detectable ground motion. Those relations then become 
equations for the limiting sensitivity of a seismic detector. It is seen 
that (7), (7’) can be used to obtain the minimum suspended mass of 
a seismic detector which is to record ground motion of a given ampli- 
tude. 

The determination of G,, Gz requires some consideration of the 
frequency spectrum of the Brownian motion. For this purpose, let the 
supposition be made that the motion is different from zero only in 
some very long time interval + which must be taken to be much 
greater than the period to which observations are extended. The force 
F in (1), which now will be regarded as due to thermal phenomena and 
which causes Brownian motion, can be represented as a Fourier 
integral 


B= (8) 


and, since thermal forces are of random nature, no frequency can 
predominate, so that 


| F(w) | = const. (9) 


The displacement of the suspended mass of the geophone is obtained 


from (1) 


wo? — w? + 2jwuo 
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and by Parseval’s formula 
| F(w) |?-dw 
0 (wt + wot)? 
Equations (6) and (11) can be used to obtain the value of | F(w)| 


which is constant and can be taken outside the integral sign in (11). 
It follows 


(11) 


kT 
| F(w) |? = a (12) 


The effect of the frequency cut-off of the recording device is to 
change the upper limit in (11) from © to a. Thus, the mean square 
‘ Brownian motion which is recorded becomes 


4 kT dw 
(x3) 
0 (w? + wo)? 


and a simple calculation gives 


wo 0” w? 
as the factor to be used in (7). A similar procedure leads to 


Wo wo”? + 


which is to be used in conjunction with (7’). 

The relations between mass and minimum detectable ground 
amplitude established by (7) and (7’) are of fundamental importance 
for seismic detectors; they apply quite generally regardless of the 
specific design adopted. But the limiting sensitivity which can be de- 
rived from those relations is not necessarily attained; under certain 
circumstances, other phenomena, closely allied to Brownian motion, 
may lower the useful sensitivity of a seismic detector. This is often the 
case for an electromagnetic geophone which will be discussed in the 
next section. 

For the purpose of numerical illustration, suppose f=fo, T= 300° K, 
and G=1.oas an approximation; employing the c.g.s. system of units, 


MA? > 16.7 X 10-"4/f,2. 


The amplitude of ground unrest is generally considered to be of the 
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order of 10~-* cm, so that a seismic detector capable of reaching that 
level must have a suspended mass 


M = 1670/f)?. 


An instrument with a period of 1 sec. should have a mass of at least 
1670 g. A to sec. period seismograph would require a mass of 167,000 g 
for the same sensitivity. Fortunately, such high sensitivity is not 
generally desired in a long period instrument; a mass of 10,000 g will 
give a limiting sensitivity of 4X10~* cm for a 10 sec. period. 


THE LIMITING SENSITIVITY AND THE POWER OUTPUT 
OF AN ELECTROMAGNETIC GEOPHONE 


The moving coil electromagnetic geophone consists of a length of . 
conductor capable of motion at right angles to the lines of force of a 
magnetic field. The ends of the conductor are connected to an am- 
plifier and a recorder; an external resistance is shunted across the 
moving conductor to provide electromagnetic damping. Let the re- 
sistance of the geophone be denoted 7, the externa] resistance R, the 
magnetic field H, the length of the conductor /. Since the geophone 
is critically damped 


PH?/(r + R) = 20M (15) 


and the equation of motion is the same as (1). 
The voltage output of such a geophone is 


R 
r+R 
and it follows from (5) and (15) that the mean square voltage can 
be written 
R? w® 
V2 = wMA? . (16) 
r+R (w? + wo?)? 

There are now two fundamental factors which limit the sensitivity 
of the geophone. One is the Brownian motion of the suspended con- 
ductor which was treated in Section 1. The other factor is the spon- 
taneous thermal voltage fluctuation across the resistances r and R in 
parallel, given by 


2 rR 
Vy? = — kTon (17) 


T r+R 
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where w;/27 is the cut-off frequency of the amplifier. The total voltage 
fluctuation is the sum of (17) and of the fluctuation produced by 
Brownian motion. As in Section 1, it follows 


2akT w + wo? \ 2 
MA?=2 ( ) ‘G3 (18) 
w? w? 
with 
rT 
G3 = + Go. (19) 


The first term on the right-hand side of (19) is the contribution of the 
thermal fluctuation voltage in the resistor, the second term is due to 
Brownian motion. 

Equations (18) and (19) can also be derived from the fundamental 
theorem on thermal voltage fluctuations in an arbitrary circuit. This 
theorem can be stated 


2 
Vy? = — rf p(w) -dw 
us 0 


where p(w) is the resistance component of the impedance of the input 
circuit of the amplifier which now includes the motional impedance 
of the geophone. 

Since G; is always greater than Ge, an electromagnetic geophone 
has a higher noise level than a mechanical seismograph. This is very 
pronounced when the cut-off frequency is much greater than the geo- 
phone frequency. From this point of view, optical registration is to be 
preferred whenever extreme sensitivity is desired at low frequencies. 
Another method to accomplish the same result is to make Rr. 

Substituting 


4 w = wo r=R T = 300°K 


and employing the c.g.s. system of units, 


MA* & 5.34 X 10°" 


if the Brownian motion can be neglected (f,>/0). 

For example, let a 5 cycle/sec. geophone be desired which is to be 
connected to an amplifier with a cut-off frequency of 50 cycles/sec. 
If the geophone is to detect ground motion amplitudes of 1o~* cm, 
the minimum suspended mass which must be used will be 210 g. 
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THE GEOPHONE AS A POWER GENERATOR 


The most practical approach to the problem of geophone design 
is through the direct experimental determination of recorder power 
requirements and computation of geophone power output. For, under 
ordinary conditions, factors other than Brownian motion and thermal 
voltage fluctuation limit the useful sensitivity of a geophone. Such 
factors are, for instance, power line and other electrical pick-up, 
microphonics in the amplifier, transformer noise. The geophone can 
then be regarded as a generator which must furnish a certain minimum 
power to the recording device, e.g., the input circuit of an amplifier. 

If the noise level be denoted «, the power to be generated by the 
geophone is ae. Now, the maximum power which a critically damped 
geophone can transmit to the recording apparatus is 


(ww? + w?)? 
because this expression represents the power, derived from ground 
motion, which is dissipated in the damping device. It follows 


2 2\ 2 
MA? 2 ae : (=). 


wow? w? 


Wo 


An amplifier of the type usually employed in geophysical explora- 
tion for the recording of reflections will operate satisfactorily with an 
input of about 5X10~ watt=5X10~° erg/sec. (this is equivalent 
approximately to 25 microvolts across a resistance of 1 megohm). 
Substituting this value of power for ae in (20), and assuming w= wo, 


MA* & 8.1 X 10°" 


Since the frequency of the usual reflection is of the order of 40 cycles/ 
sec., and the amplitude of ground unrest may be taken as 10-8 cm, 
it follows that the suspended mass of a reflection geophone should be 
greater than 12.5 g, a requirement which is easily satisfied. 


INTERPRETATION OF THE TRANSIENT BEHAVIOR 
OF THE REFLECTION SEISMOGRAPH* 


R. G. PIETYt 


ABSTRACT 


This paper discusses a method of design for the reflection seismograph based on the 
galvanometer response when the geophone is subject to a unit impulse. A graphical 
method of obtaining the actual response to an arbitrary geophone motion in terms of the 
impulsive response is given. A family of desirable types of impulsive response is obtained 
by analyzing the implications of this graphical computation when strong low and high 
- frequency interference is present. The application of this method of computation to the 
evaluation of corrections in the apparent arrival time of reflections obtained with dif- 
ferent instruments is outlined. 


It has been recognized for some years that the performance of the 
reflection seismograph depends on the overall tuning of the system. 
This arises from the fact that the motion of the ground at the surface 
is not due to reflected waves alone but contains other undesirable 
types of disturbances. The advantages obtained by tuning depend on 
the fact that the interfering ground motion, which is most frequently 
caused by wind and surface waves set up by the shot, is composed of 
either more rapid or more slowly varying oscillations than the ground 
motion caused by reflected waves. The actual reflected wavelet may 
contain only two or three oscillations so that its frequency spectrum 
is fairly broad. Hence we cannot sharply tune the instruments to it 
as in the case of sinusoidal oscillations. Nevertheless it has been found 
possible to obtain some separation of the desired waves from the un- 
desired waves by observing the average period of the oscillations in 
a reflected wavelet. The overall response of the seismograph is then 
made a maximum to sinusoidal excitation of this period. By experi- 
mental adjustment of the rate at which the response falls off on each 
side of the maximum a balance may be obtained between the degree 
of objectionable alteration in the character of the reflections by the 
tuning and the suppression of unwanted disturbance. 

The above method of design is based on the steady state behavior 
of the instruments. Most of the literature on seismograph instruments 
describes their behavior in terms of their steady state response to 
sinusoidal excitation. There are two objections to this steady state 


* Presented at the Annual Meeting, Houston, Texas, April 2, 1941. 
{ Phillips Petroleum Company, Bartlesville, Oklahoma. 
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point of view. First, the instruments are characterized by two de- 
pendent variables, the amplitude and phase as functions of the fre- 
quency of excitation, while the seismogram shows the ground motion 
as a function of another variable, the time. Second, the relation be- 
tween the amplitude and phase versus frequency and their effect on 
the record is too complicated to give more than qualitative assistance 
to the designer. This difficulty arises because Fourier Integrals are not 
easy to interpret except for the case where distortion-free reproduction 
is desired, but here distortion of the true ground motion is the essence 
of the problem. 

The purpose of this paper is to point out certain relations between 
the appearance of a seismogram and the transient behavior of the 
instruments to impulsive excitation. These relations may be used as 
the same type of guide in seismograph design that the frequency 
response furnishes in audio frequency communication design. The 
natural tool for treating transient behavior is provided by the opera- 
tional calculus. However, in many cases where laboratory facilities 
are available and the cost of experimental equipment is low, a satis- 
factory design may be developed by ‘cut and try’ methods in less 
time than by detailed computation. In this paper a number of useful 
relations which may be used to guide this type of design will be dis- 
cussed. The relations are rather elementary to appear at this late 
date, but the writer has not seen them explicitly discussed in the 
literature on the seismograph. 

The assumption is made that the behavior of a seismograph from 
geophone to galvanometer is described by a set of linear differential 
equations with constant coefficients. If a suitable member of the set 
of differential equations is equated to a time function which has zero 
value for all times less than the time zero and unit value for all times 
greater than zero, then the complete solution of the set, giving the 
galvanometer deflection, will consist of the sum of exponential func- 
tions of the time plus a constant. If A(#) is this function, then: 


A(t) = a+ t>o 
1 (1) 


A(t) =0 ix~«o 
ax OF ax may be real or complex constants. If they are complex they 
occur in conjugate pairs. The special case solutions for particular 


values of the constant coefficients do not have physical significance 
here, because the seismograph components are not measured with 


| 
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mathematical precision, and so Eq. (1) is in practice completely 
general. 

Now it can be shown that the solution of the original equation in 
terms of A(é) may be obtained for any type of time function, which is 
zero for all negative times. A convenient form for this solution is 
given by! 


i(t) = e(t)A(o) + f “Al — A)e(A)dr (2) 


where A’(¢#) means (0/d#)A(t—X). The integration parameter dis- 
appears when the limits are substituted after integration. The proof of 
Eq. (2) consists in substituting back in the original equation and veri- 
fying that the boundary conditions are satisfied. This equation is an 
analytic statement of the familiar superposition theorem. 

For a seismograph of the type ordinarily used in reflection work, 
i(¢t) may be conveniently chosen to represent the galvanometer de- 
flection, and e(¢) to represent the force acting on the geophone. With 
this choice of variables A(¢) represents the galvanometer deflection 
due to a suddenly applied unit force on the geophone at the time ¢=o. 

Fortunately, the fundamental Eq. (2) contains easily measured 
quantities. A(#) may be measured with sufficient accuracy for most 
purposes by suddenly removing a weight from the center of gravity 
of a shaking table supporting the geophone. For example, a shielded 
electromagnet may be used to release a small ball. A’(#) may be 
measured by giving the geophone a sudden change in momentum. 
In this case a small ball may be bounced on the shaking table. By 
this simple operation it is possible to obtain the functions A(é) or 
A’(t), which completely determine the entire system. Consequently, 
Eq. (2) makes a convenient starting point for the study of the transient 
behavior of the seismograph. 

In Eq. (2) the first term contains the factor A(o) which is the 
galvanometer deflection at the instant of application of a unit force. 
This term in seismograph design is usually zero. In the integral, the 
time appears as the upper limit of integration and as a parameter 
under the integral sign. This makes the interpretation more complex 
than the usual definite integral. It should be noted that at each instant 
of time ¢, A’(t—A) can only have a value different from zero if \<¢#. 
A'(t—X) will reduce the integrand to zero wherever \>#. Also, as A 
increases from zero to ¢ the argument of A’(¢—)) decreases so that the 


1 V. Bush: Operational Circuit Analysis. Wiley & Sons, New York, 1929, page 68. 
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galvanometer deflection i(¢) at a time ¢ contains a contribution from 
the ground motion for all past time. The contribution from past 
ground motion varies according to the historical weighting function, 
A'(t—X). For instance the contribution to the deflection from t=a 


to t=a+Aa at the time ¢ is 
A'(t — a)e(a)Aa (3) 


where Aa is sufficiently small. In order to visualize Eq. (2) more 
clearly, consider the seismogram in Fig. 1. At (a) we have a single 


Force on Geophone 


A> e(A) 


(a) 


A( » Transparent sheet 


(>) 


(c) 
A(t -a)e(a) 


Galvanometer Ceflection 


Fic. 1. Graphical computation of: 


The transparent sheet shown at (6) is laid on the curve giving the force on the geophone 
shown at (a). The corresponding ordinates are then multiplied and the area under the 
product evaluated as shown at (c). The curve (d) then is plotted for a series of such 
evaluations made at typical points as abcd etc. 


trace with an idealized reflection showing the actual force acting on 
the geophone. At (6) we have a transparent sheet with A’(A) drawn 


| 
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on it with A increasing to the left. Now if A’(A) is laid upon the seismo- 
gram as at (c) and we call the position of the zero axis on the sheet ¢, 
and time on the seismogram X, then obviously the curve on the trans- 
parent sheet becomes A’(¢—X) when the d axis is used as a reference. 
If at (c) the ordinates of the corresponding points of the superposed 
curves are multiplied and the product plotted, the lower curve in (c) 
is obtained. At the time ‘“’ the area under this curve is the galvanom- 
eter deflection for the time ‘.’ The deflection for any other instant of 
time is computed in a similar manner so that for a sufficient number 
of such computations as for points a, 5, c, d, the complete galvanometer 
deflection may be obtained. The curve (d) is the actual galvanome- 
ter deflection which is obtained by a seismograph, with an impulsive 
response A’(t), excited by a ground force e(#). This cinemagraphic 


A 


(a) 


(5) 


AX 


(4) 


Fic. 2. Basic forms of impulsive response A’(#). The rectangular forms at (a), (6), 
(c), (d), correspond to ideal circuits which, when properly taken to a differential limit, 
respectively reproduce the force on the geophone or take the first, second and third 
derivative of the force on the geophone. The smooth forms are physically possible 
finite equivalents for slowly varying forces on the geophone. 


process shows that when A(o) is zero, then A’(#) completely charac- 
terizes the system. If various sine waves are plotted as the assumed 
ground force, the steady state frequency and phase characteristics can 
be obtained. There are, of course, easier methods for this special case. 


! | 
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The graphical evaluation of Eq. (2) makes it possible to under- 
stand several important relations between A’(t) and e(#) which are 
essential in a well designed seismograph.’ First, it is clear that two 
reflected wavelets must be separated by an interval greater than the 
interval for which A’(t) has appreciable magnitude, in order for the 
wavelets to appear on the seismogram with an undisturbed interval 
between them. Obviously, A’(#) should occupy as short an interval as 
is consistent with other requirements, in order to give high resolving 
power. This requirement generally eliminates the use of lightly 
damped circuits or filters with a sharp cut off between the attenuated 
and pass bands, since in this case A’(¢) is unnecessarily of long dura- 
tion. 

Second, consider the case where the ground motion is slowly vary- 
ing and of long duration compared to A’(#). The result is a seismogram 
whose character is largely dependent on the number of oscillations in 
A'(t). In Fig. 2 several idealized forms for A’(¢) are shown; beside them 
there are certain typical forms which, in the limit for sufficiently slow 
ground movement, are equivalent to the rectangular forms. These 
equivalent forms for A’(#) might be observed on an oscillograph in 
connection with a study of different types of electrical networks. Con- 
sider an A’(¢) of the form shown in Fig. 2(a). In this case A’(é) will 
reduce the product function to zero everywhere except over a small 
interval AX long. Assuming A(o) =o (this will be assumed in all cases 
to follow), Eq. (2) reduces to the form, 


t 
= he(d)dX = e(t)AXh = ae(t) (4) 
t—Ad 

where a=hAd and A) is sufficiently small. This type of A’(#) merely 
reproduces the function, that is i(¢) and e(#) are proportional. If Ad is 
made large a sort of running average of e(#) is obtained for the gal- 
_vanometer deflection, i(¢). When A’(é) is modified to have the form 
shown to the right in Fig. 2(a), the running average is weighted. 
That this modified form for A’(/) gives an attenuation of the higher 
frequencies for sinusoidal ground motion, can be verified by the 
method shown in Fig. 1. It is obvious that the high frequency attenu- 
ation arises from the averaging brought about by the form of A’(é). 
2 It is implicitly assumed that the seismograph is to operate under conditions of 


strong interference. The less restricted requirements on A’(t) when only mild interfer- 
ence is present will not be discussed in this paper. 
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If A’(t) does not start at =o then i(é) will be delayed or shifted along 
the time axis. 

In Fig. 2(b) there is shown a second type A’(#). In this case another 
rectangle has been added on the negative side of the axis to the single 
rectangle of Fig. 2(a). This form of A’(#) will reduce Eq. (2) to the 
form ‘ 


i(t) = f — he(d)dd (5) 


t—2Ad 
= hA*de'(t) = be’ (t) (approximately), 


where =A?) and Ad is sufficiently small. For this type of A’(é), i(#) 
is proportional to the derivative of e(¢). Again, if AX is large or A’(t) 
is suitably modified, the result is some type of running average deriva- 
tive of e(¢) for the galvanometer deflection, i(#). It is obvious that in 
this case the high and low frequencies will be attenuated when sinu- 
soidal ground motion is present. Since at any instant the average 
ground motion for a short interval is subtracted from the motion of an 
adjacent interval, it must be true that for sufficiently slow ground mo- 
tion the galvanometer deflection becomes small. When the motion is 
very rapid the galvanometer deflection becomes small because of the 
averaging process. 

This same type analysis can be applied to an A’(#) as shown in 
Fig. 2(c) or (d). In Fig. 2(c) the A’(¢) takes a running average second 
derivative and in Fig. 2(d) a running average third derivative. Any 
actual A’(¢) which is not a generalization of the types shown may be 
considered to be built up by a suitable addition of these elementary 
types. It is possible to generalize from Fig. 2(a) backward to obtain 
A'(t) functions which integrate e(é). 

The relationships between the form of A’(t) and the response of a 
seismograph to slowly varying ground motion given above, show how 
the number of oscillations and their form affect the low frequency 
response. 

The curves at a, b, c, d in Fig. 2 taken progressively are obtained 
by a suitable addition of one half cycle to the preceding curve. Con- 
sidering the steady state frequency versus amplitude characteristic 
of the seismograph, it is evident that it is possible to obtain six deci- 
bels of attenuation per octave, as a limit, on the low frequency side, 
for each half cycle after the first in A’(¢). No general theorem relating 
to the number of cycles and the attenuation is possible when circuits 
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like the Wein bridge are to be included, but the curves of Fig. 2 show 
what can be done in the case where a steadily falling response on the 
low frequency side is required. An advantage of this purely transient 
point of view is apparent when it is considered that by directly meas- 
suring A’(¢t) any undesirable form for A’(#) may be quickly and cer- 
tainly detected and appropriate design changes made. 

Let us consider a third set of relations between A’(#) and the seis- 
mogram. If the seismograph is excited by rapidly varying ground mo- 
tion when compared to the variations in A’(#), then a point of view 
inverse to that above can be taken. In Eq. (2) and the graphical 
equivalent given in Fig. 1, e(A) and A’(#) can be interchanged, for by a 
simple change of variables Eq. (2) becomes 


i(t) = e(t)A(o) + f — A)dr. (6) 


The change in wording necessary to make all that was said in the para- 
graph above apply to this new case is elementary. The short period 
disturbances have in (6) become the historical weighting function for 
A'(t). It is apparent that in order to suppress short period random dis- 
turbances, it is necessary to keep the slope or derivatives of A’(t) 
small. In other words A’(#) must not have any sharp bends and should 
oscillate as slowly as possible. 

The reader is no doubt aware of the fact that there are bridge cir- 
cuits which will suppress a single sinusoidal wave completely. How- 
ever, if the bridge will completely suppress a wavelet of short duration 
it will be in balance at all frequencies, and therefore useless for seis- 
mograph tuning. The requirement that the derivative of A’(¢) be as 
small as possible does not give suppression for impulsive excitation of 
the form shown in Fig. 2(a). There is not much that can be done about 
this excitation since A’(#) must be such that it will have good response 
to wavelets of longer period. 

The analysis of the tuning requirements up to this point has not 
completely restricted the A’(t). The form may be varied within limits 
to make it easy to determine the time of arrival by observing the 
reflections on the seismogram. The field conditions are so variable that 
an exact analysis of the relation between A’(¢) and the probable error 
in picking the arrival time of a reflection, would probably involve as- 
sumptions that are not constant from one location to another. How- 
ever, it seems reasonable to assume that the shorter the total duration 
and the more peaked the envelope of a wavelet, the easier it will be 


/ 


TRANSIENT BEHAVIOR OF REFLECTION SEISMOGRAPH 131 


to determine its time of arrival. It is easy to see that this requirement 
means that A’(t) should be of short duration and have a peaked en- 
velope. The A’(#) shown in Fig. 2(a) meets this requirement but it 
does not adequately suppress interference. Since it is assumed here 
that the problem is to obtain reflections in the presence of quite severe 
interference, some other functional form from A’(#) must be used. For 
instance if the sharp corners of the rectangles in Fig. 2(b), (c), (d) are 
rounded, this requirement of a peaked envelope will be satisfied. 

The graphical process of Fig. 1 is quite useful to the designer in 
estimating the apparent shift in arrival time caused by changes in 
tuning. It is only necessary to make graphical computations for a few 
assumed wavelets using the directly measured A’(t). The relation 
between the times for different cycles on the resultant and on the orig- 
inal wavelet can then be evaluated. If these computations are carried 
out, it will be found, that if the original wavelet contains ‘n’ half 
cycles, that each successive derivative will add one half cycle to the 
original wavelet. (This follows analytically from Rolle’s Theorem when 
it is applied to functions of the type given by observed wavelets.) The 
added half cycle cannot be considered to have been added at any point 
in the wavelet so that in order to adjust the time correction when a 
circuit change is made it is almost essential that some equivalent of 
the graphical process of Fig. 1 be used, or that the correction be ob- 
tained by comparison of seismograms taken under identical field con- 
ditions. 

The entire discussion above makes it possible to understand the 
nature of the more usual compromises necessary to satisfy the reflec- 
tion seismograph design requirements. The requirements on A’(é) 
which will enable a seismograph to alter the apparent ground motion 
so that the arrival time of reflected wavelets may be observed in the 
presence of both short and long period interference are summarized 
as follows: 

(a) A’(¢) should oscillate several times with a short period to sup- 

press slowly varying disturbances. 

(b) A’(#) should have appreciable magnitude for only a short 
interval of time in order to give good separation between re- 
flected wavelets. 

(c) The derivatives of A’(#) should be small to suppress short 
period disturbances. 

(d) The envelope of A’(#) should be peaks to make it easy to pick 
the arrival time. 


Naturally these requirements are mutually incompatible. The best 
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balance between them depends on highly variable field conditions. 
The nature of the problem is such that it cannot be precisely stated 
for analytical treatment. However, it has been observed that some 
member of the family of curves shown in Fig. 2 (to the right of its 
corresponding rectangular equivalent) will balance the requirements 
(a), (b), (c), (d) for most field conditions. On the low frequency side 
the attenuation given by (b), (c), or (d) will increase six decibels per 
octave for each half cycle added to the reflected wavelet as it appears 
on the seismogram. This high slope makes it necessary to have only 
three or four half cycles in A’(¢#) under most conditions. If the zeros 
of A’(t) are evenly spaced at about the same interval as the zeros of a 
reflected wavelet, the derivatives of A’(#) may be made sufficiently 
small to give good high frequency attenuation. The envelopes are 
peaked so that the reflections will in turn be peaked. Actual field 
experience shows that a rather large tolerance to variations from the 
exact form of A’(t) shown in Fig. 2 is permissible. For instance an 
A’(t) which has one extra half cycle, which may be of little use because 
of its long period, will make a deficiency which in the average case is 
just noticeable. 

At this point it would be natural to proceed with a discussion of 
the type circuits which will approximate the desired forms of A’(é). 
This problem can be approached by considering the Heaviside expan- 
sion theorem of the operational calculus. It is the author’s intention 
to discuss the application of the operational calculus to this problem 
in a subsequent article. 
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CURVED PATH METHODS APPLIED TO VERTICALS 
AND TO WIDE SHOT SPREADS* 


MORTON MOTT-SMITHf 


ABSTRACT 


The ray path-wave front chart and similar methods used for vertical shots are 
briefly described, and the possibility and the advantages of applying the linear time law 
to these methods is pointed out. A short method of fitting the linear depth law to the 
observed velocity is described. The possibility of extending any method applicable to 
verticals to wide shot spreads is shown, and a method of correcting straight path dips 


for curvature in the case of verticals is pointed out. 


VERTICAL SHOTS 


It is not difficult to apply curved path methods to vertical shots, 
for which the receptors are symmetrically arranged on either side of 
the shot point. The simplest way to do this is to construct a ray path- 
wave front chart, such as recently described by Peterson,! the essen- 
tial features of which are shown in Fig. 1. 


Sind = Error varies as 
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+ Formerly with Independent Exploration Co., Houston, Texas. Now with Science 
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The path and time equations are usually given in terms of a param- 
eter, p=sin ¢/V, where ¢ is the angle between the ray and the vertical 
at any point, and V is the velocity at that point. This is constant for 
any particular ray, hence distinguishes one ray from another. Here 
it is more convenient to use for this purpose ¢o, which is the value of 
¢ at S, and will be called the starting angle. Also, since only the ter- 
minal points of the rays come into consideration, ¢ (without subscript) 
will here denote the finishing angle of the ray, such as its value at a 
reflection point, F. In the parametric equations, then, x, z, and ¢, are 
expressed in terms of ¢o and ¢. 

For a vertical shot, the middle or axial ray of those falling on the 
plane P is reflected directly back to S. The incident and reflected rays 
coincide. The starting and emergent angles are the same. The former 
can be obtained by the well-known approximation, sin ¢9= V,AT/d, 
where 6 is the receptor spread. This approximation is very accurate, 
but the error varies approximately as 6?/z?. However, if this ratio is 
large, an approximate correction can be applied. From this formula, 
the starting angle corresponding to the observed AT can be calculated, 
and thus the appropriate ray selected. Where this ray intersects the 
wave front corresponding to the observed 7, is the reflection point F, 
whose coordinates x and z can be read directly from the chart. 

The plane P that accomplishes the proper reflection must be per- 
pendicular to the ray and therefore tangent to the wave front at F. 
Consequently the dip angle 6 is equal to ¢. This angle must be calcu- 
lated from a more or less complicated equation (depending on the 
velocity law used), so that it is customary to provide a separate chart 
from which it can be read directly. 

If the linear depth law, V = Vo+ Az, is used, both the rays and the 
wave fronts are arcs of circles. The centers of the latter all lie on the 
z axis. Their positions and the radii can be calculated and the wave 
fronts put in by compass. Thus C,, is the center, and R,, is the radius 
of the wave front drawn in Fig. 1. If the centers are marked and num- 
bered on the chart, the plane can be drawn directly on the profile, 
after the manner of an ordinary dip chart, by simply making it per- 
pendicular to R,, at F. Used in this way, the chart is almost as fast as 
a straight path dip chart, and is read to the same accuracy. 

One advantage of this ray-wave chart is that if constructed for a 
particular pair of velocity constants, Vo and A, it can be used for any 
other pair (belonging to the same velocity law) by merely changing the 
x, 2, and ¢ scales. Or a new chart can be calculated by merely multiply- 
ing x, z, and ¢, by the appropriate factors. The charts are also easy to 
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calculate and construct. There is in fact nothing easier in the whole 
realm of curved path operations. 

One disadvantage is that if higher accuracy is desired, and sought 
by increasing the scale, the chart becomes both wide and long and 
therefore unwieldy. This inconvenience can be avoided by changing 
the coordinates. Thus if T and AT are used in place of x and gz, the 
AT axis will be considerably shorter than the T axis. The chart be- 
comes long and narrow and easier to handle.? Other variations are 
possible. 

Such charts are much used, but for the most part they are based on 
the linear depth law—which the author has called the mathematicians 
delight. However, they all can without much additional computation 
be constructed for the linear time law, V = Vo+At,, where #, is the time 
along a vertical ray. This law, as pointed out by the author in a 
previous paper,? fits a large range of observed velocities much better 
than the linear depth law does, and moreover can be closely fitted to 
any one of them by the method of least squares. This cannot be done 
for the linear depth law, because in the average velocity equations the 
slope constant A appears both algebraically and as part of a logarithm 
or of an exponent. Consequently these equations can only be solved 
for A by a method of successive approximation. About the best that 
can be done in fitting the mathematical curve to the observed average 
velocity curve is to make it intersect the latter at two chosen points. 
Thus only two points of the observed velocity are actually used, the 
others merely aiding in the preliminary smoothing of the curve. For 
linear time, on the other hand, all of the observed points can be used. 
A preliminary smoothing is not necessary or even desirable, because 
the least squaring process does this automatically. 

However, for linear depth, it is possible to avoid the difficulties 
with A by constructing a curve from which A can be obtained for all 
velocities that can be adequately represented by this law. Such a 
curve can be calculated directly, that is, without successive approxi- 
mations. Or, a simple algebraic formula can be set up from which A 
can be obtained approximately, but with sufficient accuracy. Once A 
is found, Vo can be calculated by a direct formula. 

In a discussion of the linear time law by Houston,’ he remarked: 


2 The manner of calculating such a chart is given by J. J. Jakosky, “Exploration 
Geophysics,”’ p. 490, Los Angeles, 1940. 

3M. Mott-Smith: Gropnysics, Vol. 4, No. 1. Jan. 1939. 

4 Charles E. Houston, “Seismic Paths, Assuming a Parabolic Increase of Velocity 
with Depth,”’ Geopuysics: Vol. 4, No. 4, Oct. 1939. 
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“As Muskat pointed out, because of the transcendental dependence 
of x and ¢ on the parameter, one cannot obtain x(¢) explicitly. The 
example is of interest however in that the equation of the wave paths 
can be obtained (in parametric form), the time-distance curve of the 
direct wave can be computed, and the corresponding reflection time- 
distance curves can be gotten by graphical methods.” This means that 
a ray-wave chart, and the other charts used for verticals, can be cal- 
culated for linear time by direct means, i.e., without approximations. 
The x equation is to be sure more complicated than for linear depth, 
but the time equation is much simpler. It contains no logarithms or 
exponentials. 

On the other hand, the linear depth law is not entirely free of trans- 
cendental relations. In fact, the time equation is transcendental, and 
this causes trouble right at the start (where trouble is most discourag- 
ing), namely, in fitting the mathematical to the observed velocity 
curve. It is also the chief trouble maker all the way through. The path 
is indeed simple, being a circular arc, while for linear time the path is 
a cycloid. 

These two circumstances somewhat offset each other, so that on 
the whole it is not much more work to compute a chart for linear time 
than it is for linear depth. 

There is in fact a curious analogy between these two velocity laws 
and those that govern the circular and cycloidal pendulums. For the 
former, the path is simple and circular, but the time of oscillation i isa 
complicated matter, being given by 


T = + sin? 6/22 + 3? sint 0/2?- 4? 
+ 32-5? sin® 6/2?-42-62 + ---), 


where L is the length of the pendulum, and @ is the amplitude angle. 
On the other hand, for the cycloidal pendulum the path is complicated, 
but the time of oscillation is very simple, being merely T=4rv/1/g, 
where r is the radius of the rolling circle. 

It may be interesting also to point out that the linear time law is 
identical with the law of falling bodies, while the linear depth law is 
the alternative which Galileo proposed but rejected, because of its 
supposed self-contradictory nature. Galileo was wrong about the self- 
contradiction, but right in another way. His difficulty arose when he 
attempted to find out what the average velocity over a given distance 
would be—the only thing he could measure. Since this velocity is 
V=Az/In(Az/Vo—1), it is no wonder that he did not discover it in 
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the days before the calculus was invented. His reasoning led him to 
conclude that whatever distance the body had already fallen from 
rest, it would have used up all of the available time, so that the re- 
mainder of the journey would have to be accomplished in just no time 
—which seemed to him absurd. However, this reasoning is not as 
wrong as it seems. Had Galileo pursued it a little further—backwards 
—by making the initial distance traversed smaller and smaller, he 
would have seen that a body under this law cannot start from rest at 
all—and this is true. If one puts Vo=o in the equation above, the 
average velocity becomes infinite. However, there is nothing to pre- 
vent somebody giving the body a push, after which it can go on pre- 
cisely as the law requires. Mechanical motions can and do take place 
in accordance with this law, provided that the body concerned has 
acquired an initial velocity by some other means. 

The same is true for seismic waves. If there is to be a finite average 
velocity, there must be a finite initial velocity. On the other hand, the 
linear time law is self-starting. Bodies do fall from rest. For the pur- 
poses of a ray-wave chart, however, an initial velocity is necessary in 
either case, for otherwise all rays start vertically downward at the 
shot point, and there is no observable quantity there to distinguish 
one ray from another, or to determine where it will come up. In fact, 
no physical reason can be given why it should ever come up, or deviate 
in the slightest (except by reflection) from the vertically downward 
path on which it started. 

For linear time the paths are; as mentioned, cycloids. The wave 
fronts are oval in shape. A ray-wave chart, as indicated by Houston, 
can be directly calculated. The reflecting planes are still tangent to the 
wave fronts at the reflection points, but since the wave fronts are no 
longer circular, the simple graphical method of determining dips, de- 
scribed for linear depth, can no longer be applied—unless one can 
draw eye tangents with passable accuracy. However, the oval wave 
front can be replaced by one or more circular arcs, with sufficient 
accuracy (except for very small 7), and the centers being given, the 
graphical method can then be applied to this approximate linear time 
chart. 

The other types of charts used for verticals can also be directly 
computed for the linear time law. 


WIDE SHOT SPREADS 


The charts for verticals can be used for end shots, and for short 
spread shots, but the shot length cannot much exceed 1000 feet with- 
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out the errors becoming excessive. For greater shot lengths, other 
methods or charts must be used. 

The problem of the wide shot spread is much more complicated 
than that of the vertical shot. Its general fetures are shown in Fig. 2, 
where X is the shot length, ¢ and ¢, are the times along the incident and 
reflected rays respectively, po and p are the starting and finishing angles 
of the reflected ray. 
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Integrable equations for #, ¢,, x, and z, have been given by Ruther- 
ford’ for a velocity that varies as the square root of the depth, which 
however, is not quite the same as the linear time law discussed here. 
It is not difficult to set up such equations, but they are altogether too 
complicated for practical use. The real problem is to devise a simple 
approximate method or chart by which results can be obtained quickly 
and with sufficient accuracy. Also the method of computing the chart 
must not be too long. The method devised by the author can only be 
given here in outline. 

First, the emergent angle of the reflected ray is computed from 
sin pp=V,AT/b. A set of rays could then be calculated for specified 
values of sin po or of AT, but since they would be identical with those 
drawn on the vertical chart, it is only necessary to trace the latter on 


5H. M. Rutherford, “Reflection Dip-Shooting Methods in Seismic Prospecting,” 
Transactions of the American Geophysical Union, Seventeenth Annual Meeting, 1936. 
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a transparent sheet. This sheet then becomes a movable ray chart, 
which in use is placed on the profile in such a way that the rays diverge 
from G. 

Second, a separate chart of constant T lines is made for a given X. 
This can be done by assuming values of ¢ and #, such that their sum is 
T, and finding the point of intersection of the corresponding wave 
fronts, which is F. The locus of these points is the constant T line 
sought. For linear depth, the equation of this line can be obtained in 
terms of x and 2; but for linear time, each point must be computed by 
a method of successive approximation. 

This is a very long process, but it can be avoided (except as an 
occasional test may be required), not only for linear depth but for any 
other usable velocity law, in the following mannet. For straight path, 
as is well known, the constant T line for a given X is an ellipse whose 
major axis lies in the shot line. For curved path, the constant T line 
is a complex oval which, however, closely approaches an ellipse. The 
major axis is still horizontal, but lies below the ground level, and the 
lengths of the axes are of course different than for straight path. But 
these quantities can be calculated, and an ellipse, on the same axes 
substituted for the oval. Or, better still, the ellipse itself can be repre- 
sented by two circular arcs—a method often used by draughtsmen. 
This approximation is especially close when, as in the present case, 
the two axes differ but little in length. The radii and centers of these 
arcs can be found, and the arcs put in by compass. If the centers are 
marked and numbered on the chart, they can be used for the graphical 
determination of dips, precisely as described for the ray-wave chart. 
Finally, if the dips do not exceed 20 to 30° (the amount depending on 
circumstances) the major arc alone will suffice, in which case the 
operation of the chart becomes very simple. 

To construct a time line chart like this, it is only necessary to cal- 
culate two points on each time line, namely, the positions of F for 
6=o° and for 6=g0°, as shown in Fig. 3. These points being special 
cases, the computations are much simpler than for the general point. 
If the origin is taken at M, x99 is the semi-major axis of the oval, and 
Zg9 is its depth below the ground. The semi-minor axis is then 29—2Zs0. 
J is the point where the two arcs join. ; 

Since an approximation to an approximation has here been made, 
it might seem that the results should be doubly wrong. Tests have 
shown, however, that the major arc lies actually closer to the true T 
line than does the true ellipse. Only along the minor arc, where the 
dips are very high and therefore unusual, do the two errors add. 
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Theoretically, a different T line chart is required for each different 
shot length; but actually the T lines change so little with the shot 
length that only a few charts are required. The movable ray -chart 
takes care of nearly all the differences caused by a change of shot 
length. Thus, up to a shot length of 3000 feet, two charts in addition 
to the one for verticals generally suffice for all the intermediate values 
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of X. If only a small range of shot lengths is used, one spread chart will 
do. Finally, it is possible to specify in advance the maximum allowable 
depth error, and to determine the values of X for which charts must 
be made, and therefore the number of them that will meet the specifi- 
cation. 

Image or mirror point charts can also be made for curved path. 
They can in fact be derived from the above described reflection point 
charts by a graphical method. But for such charts, both the AT and 
the T lines change with X, so that a complete chart—containing both 
sets of lines—must be made for each required X. The number of charts 
required for a given accuracy, however, is the same as for the reflec- 
tion point charts, and the charts once made are used in exactly the 
same way and with exactly the same speed as straight path charts of 
the same sort. 

In all the methods described above, curved path is applied directly 
to the determination of dips and depths. When, however, the dips are 
mild, the differences between straight and curved path methods are 
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very small, and usually quite negligible. It frequently happens that 
nearly all the dips of a prospect are mild, only a few of them being 
steep enough to require a curved path method. In that case, it would 
be advantageous to have a simple rapid method by which the steep 
dips could be corrected for curvature, while the others are let stand. 
This would avoid the calculation and use of a more complicated 
curved path chart for the whole prospect, for the sake of only a few 
dips that may require it. For purposes of correction only, the fit of the 
mathematical to the observed velocity does not have to be so good as 
for the direct use of curved path, because the errors of misfit now affect 
only the corrections, not the main quantities. Hence the simple linear 
depth law will suffice. A simple correction formula for the dips has been 
found, which covers all velocities within the practical range with an 
error no greater at any point than 5 feet per 1000 feet, and usually 
much less. For this formula, the slope constant, A, of the linear depth 
law must be known. But this can easily be found by the simple rapid 
method already described. Vo is not needed. So far, however, this for- 
mula applies only to vertical shots. 

In concluding, the author wishes to thank the Independent Ex: 
ploration Company for permission to publish this paper. 


A NOTE ON A RATIONALIZED VELOCITY 
DEPTH EQUATION* 


NEIL R. SPARKS 


ABSTRACT 


An average velocity depth equation of the form V = Vo(1+aZ)/(1+5Z) is suggested 
and its properties relative to the generalized velocity variation within the sedimentary 
column are discussed. 


The representation of the average velocity as a function of depth 
by analytic equations has long been used in seismic surveying as a con- 
venient method of simplifying and speeding up routine computations. 
The simplicity of an equation and the accuracy with which that equa- 
tion fits the observed data are generally incompatible to a degree which 
depends upon the form of the equation and the characteristics of the 
particular set of observed data. 

A more or less universal characteristic of the seismic velocity data 
obtained from velocity determinations in deep wells is that the average 
velocity increases with depth more rapidly at shallower depths than 
at the greater depths. The velocities observed from refraction shooting 
with very long spreads (depth of penetration greater than that of the 
drill) rarely exceed 20,000 ft./sec. (occasionally 25,000 ft./sec.) within 
the sedimentary column irrespective of the depth of the column. The 
corresponding average velocities to these greater depths are, of course, 
less than this figure. 

A rationalized average velocity depth equation to fit these general 
characteristics would be one which has a specified initial velocity at 
zero depth, increases with the depth most rapidly at the shallower 
depths, and approaches but does not exceed a value of the order of 
20,000 ft./sec. at the greater depths. An equation of the form, V= Vo 
(1+aZ)/(1+6Z),! possesses these characteristics. When the depth 
(Z) is zero, the average velocity is Vo. As long as the constant a is 
greater than 5, the rate of increase in velocity with depth will be great- 


* Presented at the Annual Meeting, Houston, Texas, April 2, 1941. 

Tt Joint Geophysical Laboratory Stanolind Oil & Gas Company and Western Geo- 
physical Company. 

1 A rectangular hyperbola whose asymptotes are Z= —1/b and V=Voa/b. The in- 
stantaneous velocity corresponding to this average equation is v= Vo(1++-aZ)?/(1-+-2bZ 
+abZ?). 
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est at the shallow depths. When the depth is increased without limit 
the average velocity monotonically approaches a maximum value 


of Voa/b. 


The determination of the constants (a, 6, and Vo) to fit a set of 
observational data may be accomplished in several ways. A method 
found to be convenient and rapid applies equally well either to the 
velocity-depth or time-depth plot of the data. It consists of choosing 
three points on the observed curve (Vi, V2 and V3 or Ti, Tz and 73) 
such that the corresponding depths (Z;, Z; and Zs) are in the ratio of 
1:2:4 respectively. For this condition the three constants are: 
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THE PRODUCTION OF ELASTIC WAVES BY 
EXPLOSION PRESSURES. I. THEORY 
AND EMPIRICAL FIELD 
OBSERVATIONS* 


JOSEPH A. SHARPE 


ABSTRACT 


A solution to the problem of the wave motion produced when a pressure of arbitrary 
form is applied to the interior surface of a spherical cavity in an ideally elastic medium 
is derived. This solution is shown to be in qualitative agreement with a number of field 
observations of the effect of shot-point conditions on the characteristics of reflection 
seismograph recordings. 


INTRODUCTION 


Of the three physical processes involved in seismic exploration, 
namely, the initiation of the seismic waves, their propagation, reflec- 
tion, refraction, and dispersion, and the recording of some function of 
the motion of the surface, we possess the least satisfactory under- 
standing of the initiation process. This is in spite of the circumstance 
that it has been recognized from the very beginnings of reflection ex- 
ploration that the intensity, quality, and characteristics of seismic 
records are strongly influenced by conditions at the shot point, and 
despite the fact that a number of empirical observations concerning 
the generation of seismic waves have accumulated over a long period 
of field experience. 

These empirical observations, which have become a part of the 
lore of every field operator, are as follows:! 

1. A given amount of explosive detonated in a clay or water- 
saturated sand formation results in a greater amplitude of reflected 
motion than an equal charge detonated in a dense, rigid formation, 
such as limestone. 

2. If a hole is sprung by an initial large charge in order to form 
a sizable cavity, later small charges will result in a larger amplitude 


* Presented at the 1937 Spring meeting, Los Angeles. Manuscript received Feb. 10, 
1942. Part II will be published at a later date. 

¢ Joint Geophysical Laboratory of the Stanolind Oil and Gas Company and the 
Western Geophysical Company. 

1 See McReady. Gropuysics, V, 4, pp. 374, 379-381 (1940) and Beers, VI, 1 (1941), 
PP. 54-55. 
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of reflected motion than would be produced in the absence of spring- 
ing. 

3. The frequency spectrum of reflected motion is a function of the 
formation in which the charge is fired: shots fired in the low velocity 
zone result in very low frequency motion, compared to shots fired 
below the low velocity zone; shots fired in a rigid material, such as a 
limestone, result in a much higher frequency motion than shots made 
in, for example, a shale; in a general way, the high frequency content 
of reflected motion increases with an increase in shooting depth. 

4. The frequency spectrum of reflected motion is a function of 
charge size, a larger charge having a tendency to increase the propor- 
tion of low frequencies in the reflected motion. 

5. The amplitude of reflected motion produced by a given quan- 
tity of a high speed explosive, such as a dynamite, is considerably 
greater than that produced by a quantity of a low speed explosive, 
such as black powder, even when this quantity has been selected such 
that the maximum pressure attained is the same for both explosives.” 

That there has been continuing interest in the effects of shot point 
conditions on reflections records is evidenced by the issued patents 
which relate to means of controlling the initiation of elastic waves in 
order to secure reflection records of desired intensity or character. 
These patents fall into three categories: some aim, through treatment 
of the shooting cavity, to improve the efficiency of transmission of the 
explosive energy into the ground;? others describe means of selecting 
that formation as the site of the detonation which is the most likely 
to produce correlatable records,‘ and still others disclose arrangements 
which are intended to impart to the ground a motion whose spectrum 
is rich in the frequency components which are observed in the average 
recorded reflection.® 

As part of a long-range program of seismological research, we have 
investigated the mechanism of the initiation of elastic waves by explo- 
sives and the influence on reflection records of conditions at the shot 
point. This investigation has employed two lines of attack: the mathe- 
matical, with the objective of securing a rational picture of the 
mechanism of the initiation process which would explain as many as 


2 See Nash and Martin, Geopnysics, I, 2, pp. 239-256 (1936), for a discussion of the 
properties of explosives. 

5 Johnson, U. S. $2,028,286; Weatherby, U. S. $2,049,236; Salvatori, U. S. #2,073,- 
493 


4 Weatherby, U. S. #2,154,548; Faust, U. S. #2,156,624; Innes, U. S. #2,192,973. 
5 Scherbatskoy, U. S. $2,156,198; Green, U. S. $2,203,140. 
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possible of the empirical observations which are listed above, and 
would serve as a guide to experiment; and the experimental. Part I 
of this paper discusses the analytical investigation, and Part II, to 
appear in a subsequent issue of this JOURNAL, describes the results of 
some of the preliminary experiments which were performed. 


ANALYSIS 


The problem of the generation of elastic waves by explosion pres- 
sures is amenable to mathematical analysis only if one makes an ideal- 
ized statement of the problem, ignoring many factors which are pres- 
ent in actual shooting. I have chosen to state the problem as follows: 
Given a spherical cavity of radius a within a homogeneous, ideally elastic, 
infinite medium of density p and compressional wave velocity v; to find 
the elastic wave motion which results from application of an arbitrary 
pressure p(t) to the interior surface of the cavity. Fig. 1 illustrates the 
statement of the problem. 


AM,p 


Fic. 1. Illustration of the idealized model of the problem of the generation 
of elastic waves by explosion pressures. 


We should note the extent to which the problem as stated departs 
from the conditions commonly encountered in practice. The actual 
site of the explosion, which is a bulbous, water-filled cavity with a 
cylindrical, water-filled drill hole leading into it, has been replaced by 
an empty, spherical cavity. This substitution is probably justified: 
the cross-sectional area of the drill hole is generally small compared to 
that of the cavity; departure of the cavity from exact sphericity is not 
likely to produce any essential change in the result; neglect of the pres- 
ence of the tamping fluid is probably not serious, since experimentally 
we have found that when an explosive is fired in water, the form of the 
pressure wave is maintained with increasing distance from the source, 
and the pressure is inversely proportional to the distance. 
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By assuming a homogeneous, infinite medium, we have excluded 
from consideration reflections and refractions at elastic discontinui- 
ties near the source. This simplification does not impair the applica- 
bility of the analysis, since the only effect of the presence of these sur- 
faces of discontinuity is to add to the primary transient emanating 
from the source certain other components, delayed in time, and pos- 
sibly modified in form. 

By assuming the cavity to be hollowed out of an ideally elastic 
medium we have automatically precluded any consideration of the 
undoubtedly important effects of fracturing, crushing, and compaction 
of the material surrounding the cavity,® and also of the effect of the 
viscosity and internal resistance of the medium on the form of the 
wave. The assumption of perfect elasticity is equivalent, in part, to 
supposing that the cavity is so large, or the charges employed are so 
small, that the elastic limit of the material surrounding the cavity 
is never exceeded. 

To solve the problem which has been formally stated above we 
seek an expression for the elastic wave displacement which is a solution 
of the elastic equation of motion, vanishes at infinity, and satisfies 
the boundary condition that the stress in the medium must reduce, 
at the surface of the cavity, to the pressure applied to the interior 
surface.’ 

The elastic equation of motion is 


pd?q/dt? = (A+ 24) 


in which q is the displacement, and \ and yu are the Lamé elastic con- 
stants, uw being the rigidity. A solution of this representing compres- 
sional wave motion is 


q= V4, (1) 
provided the potential ©® satisfies the wave equation 
= v? 774, (2) 


6 McReady, loc. cit., discusses the results of these physical changes in the material 
surrounding the cavity. 

7 Analysis substantially equivalent to that which follows was developed for the 
discussion of a seismological problem by Kawasumi and Yosiyama, Earthquake Re- 
search Institute of Tokyo Bulletin, v. 13, pp. 496-503 (1935), and was subsequently 
utilized by Sezawa, ibid., pp. 741-749, and Sezawa and Kanai, ibid., v.14, pp. 11-17 
(1936). However, due to the limited availability of this publication, the analysis will 
be outlined here. 
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where v= {(A+ 2u)/p}"/? is the velocity of propagation of compressional 
waves. 

Polar coordinates are appropriate to this problem, and since the 
pressure is to be applied uniformly over the interior surface of the 
cavity, all motions and stresses will be independent of the angular 
variables. In view of this, (2), (1), and the boundary condition can 
be written 

0?(rb) /dt? = (3) 


u=0d6/dr, v=0, w=o; (4) 
and 
— [(A + 2u)du/dr + = p(t). (5) 


_ The only non-vanishing component of the displacement in the elastic 
wave is u, the radial component. The left side of Eq. 5 is the radial 
stress in the medium, which must become exactly equal and opposite 
to the pressure applied to the interior surface of the cavity when the 
distance is taken equal to the radius of the cavity. 

At this point it is convenient to assume \=uy.° Rewriting Eq. 5 
in accordance with this assumption, and also representing the applied 
pressure by a Fourier integral, the boundary condition becomes 


— pv?[du/dr + = 1/20 f f dy. (5a) 


A solution of the wave equation (3) which represents a wave 
diverging from the surface of the spherical cavity is 


= (1/r)eninr (6) 


in which r=¢—(r—a)/v. The right hand side of Eq. 6 multiplied by 
any functions independent of 7 and £ is still a solution: consequently 


= 


is a solution. Likewise, a sum of such solutions is still a solution, so 
that 


= ff dn dy (7) 


8 This assumption, which is substantially correct for the material below the low 
velocity zone, gives a considerable simplification of the analysis which follows, without 
modifying the final result in any essential way. 


> 
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is a solution. Forming u=0/dr and du/dr from Eq. 7, and substitut- 
ing the resulting expressions into the boundary condition (5a) gives 


— pv?/2ma f A(n) {4/302 — 4in/3av — n®/v?} dy 


This is satisfied identically provided the arbitrary function A(m) has 
the value 


A(n) = (a/p){n® + 4inv/3a — 
Substituting this expression for A(m) in Eq. 7 we get 


The potential of Eq. 8 satisfies both the wave equation and the boun- 
dary condition, and is a formal solution to the problem of the wave mo- 
tion produced when a pressure of arbitrary form is applied to the i in- 
terior surface of a spherical cavity. 

It will be desirable to secure an expression for the wave motion 
produced by the application of a pressure of unit function form to the 
interior surface of the cavity, because from this expression the wave 
motion produced by any other form of pressure application can be 
calculated conveniently by the use of Duhamel’s integral. However, 
we cannot proceed immediately to let p(#) have unit function form in 
Eq- 8, since p(#) can have a Fourier integral representation, as in 
Eq. 5a, only if the integral /*., | p(t)| det exists. Consequently, we shall 
assume, instead of unit function form, that the applied pressure is 
described by 


p(t) = for t20 


= for i<o 


in which fp is a constant representing the initial, and highest, pressure 
attained, and a is a positive constant. 

Using this form of pressure in Eq. 8, carrying out the y integration, 
and factoring the denominator, we have 


(9) 
amprd — [nia] [n+20(i+v/2)/3a] [n+20(i—v/2)/3a] 
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This integral can be evaluated by the theorem of residues, using for 
the path of integration the contour shown in Fig. 2. The result is 


Fic. 2. Contour used in evaluation of the integral of Eq. 9, showing 
the location of the poles of the integrand. 


apo/pr — —wr/V2 


sin wr + cos wr} | for r2Z0, (10) 


=o for r<o, 


in which we have introduced the abbreviation w= v/ 3a. 

We now permit a to approach zero,® i.e., we permit the applied 
pressure to take on unit function form, namely p(t)=o for t20, 
p(t) =0 for t<o. Then Eq. 10 becomes 


2a po 
3pw*r 
=o for tr<o. 


[— 1 + sin (wr + tan-! ~/2)] forr = o, 
II 


From this we get for the displacement in the medium correspond- 
ing to the application of unit function pressure to the interior surface 
of the cavity 


(=) sin (wr+tan—! v/2) 
Or 
(12) 


a = 
+Vi( sin or | for 720, and u=o forr<o, 


® The justification of this step is the verification, by substitution, that Eq. 11 
satisfies the differential equation (3) and the boundary condition (5). 


plane 
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in which we have reintroduced the definitions of w and p into the coefhi- 
cient outside the bracket. 

It is to be noted in Eq. 12 that the displacement is made up of 
three terms, one of which is non-oscillatory and two of which are 
oscillatory. The non-oscillatory term, which represents the eventual 
static displacement in the medium, and the first of the oscillatory 
terms, depend on distance as (a/r)?, whereas the second oscillatory 
term depends on distance as a/r. Consequently, the displacement at 
distances more than a few times the radius of the cavity is represented 
quite adequately by the simple expression 

“= -e-erl’2 sin wr for r2o 
2\/2ur (13) 


=o for <o. 


Thus, the theory which we have outlined predicts that the dis- 
placement in the wave which results from application of a unit func- 
tion form of pressure to the interior of a spherical cavity in an infinite, 
homogeneous, ideally elastic medium is a damped sinusoid, with an 
amplitude directly proportional to the pressure and the area of the 
cavity, and inversely proportional to the rigidity of the medium. Al- 
though the frequency f= w/24=(1/2/3m)(v/a) is directly proportional 
to the velocity of wave propagation in the medium and inversely pro- 
portional to the radius of the cavity, the damping is so high,!° 
(h=1/\/3, €=9:1), that the motion is of the nature of a pulse with a 
duration At= (32/2/2)(a/2). 

I have used Duhamel’s integral in the form 


U(r) = “ = nan, 


in which u(r) is the motion resulting from application of unit function 
pressure, and U(r) is the motion resulting from application of arbi- 
trary pressure p(t), to compute the motion corresponding to the 
several forms of pressure application which are shown in Fig. 3. In_ 
each case the form of the applied pressure as a function of time is 
shown as a dotted curve, and the displacement of the medium in the 
resulting wave motion as a solid curve. The maximum pressure at-'‘ 


10 When the analysis is carried out without making the assumption \=y, the damp- 
ing factor turns out to be h=[u/(A+2,)}!?=[(14—¢)/(1—«)]"”, in which o is Poisson’s 
ratio. For the material of the low velocity zone o may be as large as 0.4, leading to a 
damping factor about two-thirds as large as that above, which is still high. 
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tained is identical in all cases. In Fig. 3-a, the wave motion due to unit 
function pressure is shown. Fig. 3-b shows the response to pressure of 
the form p(t) = po(1—e~**) for 20, p(t) =o for t<o, where it has been 


5 to 
[\\ 
3 > 6 


NL \ 
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0 Ne = = 0 
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Fic. 3.-Displacement in the elastic wave motion (solid lines) produced by the appli- 
cation of various forms of pressure (dotted lines) to the interior surface of a spherical 
cavity in an elastic medium. Pressures and displacements plotted as functions of wt 
and wr, respectively. In all cases p(t)=0 for t<o. a.) p(t)= po. b.) p(t) = po(1—e7 v2"). 
c.) p(t) = €.) p(t) = £.) p(t) =Npole*!? 
—e-v2t), in which NV normalizes the maximum pressure to fo. 
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assumed that a=1/2w. Fig. 3-c is similar, but with a slower pressure 
build-up corresponding to a=w/+/2. In 3-d, the applied pressure is 
p(t) = poe~** for t20, p(#)=0 for t<o, in which a is assumed to be 
=w//2. Fig. 3-e is similar, but with a=+/2. Fig. 3-f shows the re- 
sponse to a pressure represented by p(#)=Npo(e-*'—e-**) for t20, 
p(t) =0 for t<o, in which WN is a normalizing factor, and a has been 
taken as a=w/./2 and B as B=+/2w. The pressure assumed in this 
case probably corresponds quite closely to that which exists in prac- 


tice. 
DISCUSSION 


We can now examine the extent to which the predictions of the 
theory which has been developed agree qualitatively with the em- 
pirical observations which were listed in the Introduction. 

Eq. 13, or equations derived from it for other choices of the form 
of pressure application, indicate that the amplitude of the displace- 
ment in the elastic wave will be proportional to a*fo/y. Consequently, 
the displacement produced by shooting in a formation of low rigidity, 
such as a clay or shale, should be greater than that when shooting in 
a formation of great rigidity, such as a limestone, which agrees with 
observation. 

It is observed that the preliminary springing of a hole with a large 
charge results in larger amplitudes of reflected motion being obtained 
from subsequent small charges than if the cavity had not thus been 
enlarged. Since the displacement is proportional to a*o, and the 
pressure wave is propagated in the tamping fluid with its amplitude 
decreasing inversely as the distance from the charge, we see that the 
displacement in the medium should be directly proportional to the 
radius of the shooting cavity," in agreement with this observation. 

Comparing Fig. 3-b and 3-c which represent displacements result- 
ing from the application of pressures which gradually increase to a 
maximum, to Fig. 3-a, representing the displacement resulting from 
a sudden application of pressure, we see that a powder producing a 
slow application of pressure should result in a smaller amplitude of 
wave motion than an explosive producing a rapid application of pres- 
sure, which is in agreement with field observation. 

According to Eq. 13, an important component in the wave motion 
resulting from the application of an explosive pressure to the interior 
surface of a spherical cavity in an elastic medium will have the fre- 


11 Edge and Laby, Imperial Geophysical Experimental Survey, p. 205, arrive at the 
same conclusion through use of an argument similar to the one which I use to explain 
the effect of charge size on the frequency of the motion. 
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quency f=w/2r=(v/2/3m)(v/a), in addition to the components de- 
pending on the particular form of pressure application. This predicts 
that shots fired in a formation of high rigidity should produce wave 
motion which is richer in high frequencies than shots fired in a less 
rigid medium, since not only will the velocity be higher in the more 
rigid medium, but also the shooting cavity will probably be smaller. 

Consider what occurs when a charge is fired at the bottom of a 
fresh drill-hole. A wave in which the stress greatly exceeds the strength 
of the material diverges from the site of the explosion, crushing the 
material as it goes. The stress in the wave rapidly attenuates, however, 
at a rate which probably is greater than inverse second power, due to 
the loss of energy by crushing. At some distance from the source which 
is a function of the strength of the material and the size of the charge 
which was detonated, the stress will be exactly equal to the strength 
of the material. At greater distances the material will be ‘“‘perfectly 
elastic,” except for viscosity and internal friction. This critical dis- 
tance, which can be called the radius of the equivalent cavity, is to be 
associated with the radius a of the cavity in an ideally elastic medium 
which enters into the theory which has been developed. We see, then, 
that a small charge, with which is associated a small equivalent cavity, 
should result in wave motion which is richer in high frequency com- 
ponents than a large charge, with which is associated an equivalent 
cavity of larger radius, in accordance with empirical observations of 
the effect of charge size on the spectrum of reflected motion. 


PROBLEMS IN TEMPERATURE CONTROL 
OF GRAVIMETERS* 


DAYTON H. CLEWELL{ 


ABSTRACT 


The theoretical considerations entering into the design of thermostatically con- 
trolled ovens for gravimeters are presented. The author points out that to obtain the 
most uniform temperature distribution without the use of excessive amounts of mate- 
rial, plural thermostat control points on the heater element are desirable. To reduce 
long-time fluctuations it is recommended that dual thermostating be employed in which 
an inner controlled-shell is surrounded by an outer thermostated layer held at an ap- 
pee constant temperature several degrees lower than the temperature of the 
inner shell. 


INTRODUCTION 
The fundamental provision of any gravimeter is a mechanism to 
translate gravity changes into increments of some directly observable 
parameter, 6. However, this parameter will be not only a function 
of gravity, g, but will also be dependent on the mass of the gravity 
responsive element, the elastic modulus of one or more springs, and 
various dimensions of the gravimeter. In general we write 


6= f(g, m, k, x”, 2"). (1) 


From this equation the gravity sensitivity of the gravimeter may 
be defined as S in the equation 


dg 
(2) 
§ 
where d@ is the increment in the observable parameter resulting from 
the gravity change dg/g. A number of the independent variables in 
Eq. 1 are likely to be responsive to a change in temperature dT, 
and a corresponding change in @ will be given by 


dT. (3) 


where >) is the overall temperature coefficient of the instrument. Con- 
temporary gravimeters have a probable error of the order of 0.01 to 
0.03 milligals corresponding to significant values of 

d 

= > 2X 


* Presented at the Annual Meeting, Houston, Texas, April 2, 1941. 
t Magnolia Petroleum Co., Dallas, Texas. 
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To maintain this order of accuracy it is immediately evident from 
Eqs. 2 and 3 that ), dT S 2X10-%. 

If it be assumed that the variables k, x, y and z are temperature 
responsive and enter Eq. 1 as factors, then 


=o+ pit (4) 


if we define o as the temperature coefficient of the elastic modulus and 
6, 8 and y as the expansion coefficients of x, y and z respectively. For 
most materials used in gravimeter construction these temperature 
coefficients are of the order of 107°. It is generally possible to arrange 
for some of the exponents #, q and r, and also some of the temperature 
coefficients to be negative so that >, itself is of the order of 10~*. This 
value of ew requires, then, that dT be of the order of 2X10~*°C. 

It is of course possible to so choose the exponents Pf, g, and r and 
the temperature coefficients ¢, 6, 8 and y in such a way that >, =o and 
one might therefore suppose that dT could be very large; i.e., no tem- 
perature control necessary. Although it is possible to adjust >> to 
zero it is still necessary to guarantee that all elements of the instru- 
ment undergo simultaneous changes in temperature. If, for example, 
the element x is located at one end of the instrument apart from the 
other critical elements, it is easily possible for x to undergo a tempera- 
ture change that may be greater or less than the temperature change of 
the other end, unless means have been provided to insure a constant 
distribution of temperature throughout the gravimeter. The usual 
practice is to use some amount of temperature compensation but the 
writer has found that even for instruments whose largest dimension 
is of the order of eight inches, nothing is gained in constancy by reduc- 
ing >, to much less than 2 or 3 X10~*. To use lower values of ),, so much 
material of high thermal conductivity must be used in maintaining 
constant temperature distribution that it seems more practical to use 
temperature control of d7, at the most, not greater than 107’, along 
with temperature compensation. A slowly changing temperature can 
be easily tolerated since it becomes a part of the instrument drift de- 
termined by repeated base-station observations. Depending on the 
degree of temperature compensation it is sufficient, then, to limit the 
temperature changes to no more than 0.002 to 0.006°C per hour. 


THERMOSTAT THEORY 
To consider briefly some of the fundamental characteristics of tem- 
perature control, assume a simple thermostated element J in Fig. 1. 
whose conductivity is k, density is p and heat capacity c. For simplic- 


a 
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ity, let J be a round rod of length / and of unit cross-sectional area, 
surrounded, except for one end, by a thermal insulating material whose 
conductivity is very low compared to k. Heat is either subtracted or 


4 
Xo 


Fic. 1. Single thermostated element. 


added to J at its uncovered end. This heat is under the control of a 
thermostat located a distance x9 from the end according to the equa- 
tion 

where x is measured along the axis of J, and Q and T are heat flow and 
temperature respectively at x at the time ¢. The symbol s denotes the 
sensitivity of the thermostat. Eq. 5 states that the thermostat tends 
to regulate the temperature to the value T=o. 

Since no heat is generated or lost except at the boundary, the equa- 
tion of continuity is applicable in the form 


A’T = — —= (6) 


and Eq. 5 is the important boundary condition to consider. 
A general wave solution of (6) is 


| 


where the frequency n= abk/pc and the velocity v= 2bk/pc, —1). 
The constants a and b are the real and imaginary parts of the complex 
constant z=a+7b which was introduced on separation of T into two 
factors involving x and ¢ separately. 

The first term in the brackets represents waves (one for each value 
of z) of frequencies m, and velocities » travelling in the positive 
x-direction, while the second term represents similar waves travelling 


oT 
ot 
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in the opposite direction. At any given time the amplitudes of the 
waves decrease in their direction of travel showing that they are space 
attenuated. However, at any given value of x the amplitudes may be 
decreasing or increasing with time, as the a’s are either greater or less 
than the 6’s. If we assume x» to be small compared to / the backward 
travelling waves reflected from the right hand end of J can be neg- 
lected because of the space attenuation factor. Accordingly, we may, 
for our simple case set all the B’s equal to zero. 

The boundary condition, Eq. 5, is applied by equating the coeffi- 
cient of each term of the series 


z=0 ax 


= kA ab(ia = b)e~ kt/peg—2anit 
0 b=0 


to each corresponding term of the series 


ST = A (a?—b*) kt/pcg—2anit (9) 


a=0 b=0 
with the result that 
seto(ia—b) = R(ia — (10) 


from which 
a 


tan 


a 


sin axo 


Since neither a@ nor b may be negative, s cannot be less than 
mk/2xo. If the sensitivity is less than this critical value, Eq. 7 degen- 
erates into the usual form of 


T = Doe’ *lee(p cos ax + q sin ax) (11) 


a=0 


and the solution is no longer of an oscillatory nature. 

Going to higher thermostat sensitivities we find that for rk/2xo 
<s<5k/2x there is but one set of values for a and 4, and 6 is always 
less than a. This corresponds to a single time-damped wave and after 


(8) 
b= — 
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a short time T reaches a steady state value. In fact b<a until the sen- 
sitivity reaches the critical value of s=35k/x»; a condition when b=a, 
and persistent oscillations occur because of the absence of time damp- 
ing. However, as long as s>52k/2xo there are two sets of values for 
a and }, the second set corresponding to a time-damped wave of 
higher frequency than the fundamental. When the fundamental wave 
becomes self-sustaining (s=35k/xo) this second wave is still time 
damped, and it becomes self sustaining only when the sensitivity is 
raised to s=141k/x9. Now, however, 6>a for the fundamental and 
the amplitude of the fundamental increases with time, limited only 
by the physical capacity of the heat source. The physical limitations 
of the heat source cause the average sensitivity to drop to the critical 
value of s=35k/xo, the harmonic wave damps to zero and we are left 
with only a self-sustaining fundamental; i.e., b=a. Thus the final 
expression for the temperature, after sufficient elapse of time to re- 
move transients, is found to be either a steady state value for s 
<35k/xo or an oscillatory solution whose amplitude is a function of 
x only, if s235k/x9. Whence 


T = (12) 


Since 6 must equal a for a steady oscillatory solution; Eq. 10 re- 
quires b=a=37/4%9 from which the frequency of the wave in Eq. 12 
is found to be 


1.75k 


(13) 


n= 


Another important feature of Eq. 12 is the space attenuation fac- 
tor. The ratio by which the amplitude of a temperature wave of fre- 
quency 1 is decreased as it passes through a material of thickness x is 


T; 
To 


= (14) 


a well known relation (1).* 7 
Eq. 8 states that the amplitude of the heat flow at «=o is Aaky/2 
and since this must equal Q,,, the maximum output of the heat source, 
in the case of self-sustained oscillations is 
Om, 
aky/2 /2mpckn 


* Refer to Bibliography at end of article. 
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showing that the amplitude, A, of the temperature oscillations de- 
creases with increasing frequency; a very desirable characteristic. 

This simple thermostated element has been sufficient to demon- 
strate the fundamental nature of thermostat action, especially for the 
case of high-sensitivity control where a steady state solution is impos- 
sible unless, by placing the thermostat in close thermal contact with 
the heat source, one can make x» so small that s<35k/ao. These char- 
acteristic features have been demonstrated both experimentally and 
analytically by L. B. Turner (2) in England. 

In the construction of any constant temperature chamber it 
should be realized that absolute constancy is impossible because for 
every change in the ambient temperature there is a corresponding 
small but finite change within the chamber. If the source of heat for 
the chamber is electrical, fluctuations in the supply voltage will also 
produce changes within the chamber. Let T be the inside temperature 
of the oven, ¢ the ambient temperature and V the voltage applied to 
the heating coils. The heat flow Q through the oven wall per unit area 


will be : 
Q=C(T 4), (16) 


where C is the conductivity of the oven wall. Thus d0=C(dT—dg@). 
Variations in Q are a result of thermostat action as well as variations 
in supply voltage so that 


dQ = ——dT + —dV. I 
(17) 
. Eq. 5 shows that 0Q/0T=-—s and since the heat generated by a 
resistance is proportional to the square of the voltage, 00/dV = 2Q0/V. 
Making the substitutions we have finally 


aQ aQ 


= + T - ——| 8) 
(1 


As long as s is large compared to C the variations inside the oven 
are of course very small, but s would need to be infinite for dT=o. 
This argument is based on the assumption that the thermostat has 
continuous control of the heat supply, Eq. 5. The argument is also 
true for relay-type thermostats which have on-off control. With relay- 
type control the thermostat is characterized by a certain temperature 
backlash (Ty—T 1) through which the temperature must oscillate. 
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If Qm is the available supply of heat the thermostat cycle will be sym- 
metrical when Q/Q,=%, and the amount by which the temperature 
T overshoots Ty is equal to the amount it undershoots 7. If, how- 
ever, T—@ increases so that the heat must be on for a considerably 
longer time than it is off, the amount of undershoot will exceed the 
amount of overshoot and the mean operating temperature of the ther- 
mostat drops below the symmetrical value of (Ty—Tz)/2. 

It has been experimentally demonstrated by Turner (3) that con- 
tinuous thermostat control is practically indistinguishable from relay- 
type control, when dealing with the very high thermostat sensitivities 
that would put the former in the oscillatory region of operation; even 


a 
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Fic. 2. Shielding effect of a cylindrical shell. 


the respective frequencies of oscillation are equal. The Turner refer- 
ence is suggested as an excellent account of the details of a precision 
isothermal chamber for laboratory operation. 


TEMPERATURE DISTRIBUTION 


In order that temperature compensation may be employed in the 
construction of a gravimeter it is essential that the temperature dis- 
tribution throughout the thermal housing be maintained as uniform 
as possible. If an ambient temperature gradient exists there will be 
a corresponding temperature gradient within the housing. The prob- 
lem is to make the ratio of the latter to the former as small as pos- 
sible. 

Fig. 2 shows a Jong cylindrical shell of high conductivity hi, 
immersed in a medium of low conductivity ke, in which medium there 
exists a temperature gradient 0¢/dx perpendicular to the axis of the 
cylinder. La Place’s equation is applicable since no heat is generated 


( 
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or lost and we are interested only in the steady-state solution. This 
problem is identical with the problem of magnetic shielding where 
k,/k2 is analogous to the relative permeability. The shielding factor of 
a single shell as represented in Fig. 2 turns out (4) to be 


4 
Of D 

‘Ox 


where d is the thickness of the shell and D>d is the mean diameter 
of the cylindrical shell. The ratio k,/ke may be as high as 10+. 

As might be expected from experience with magnetic shielding, a 
more efficient distribution of temperature with a given weight of high 
conductivity metal is obtained by dividing the metal into a series of 
concentric shells. Wills (5) has shown that the most efficient use of 
material is produced if the ratio of the thickness of the shell (including 
the shells of low conductivity as well as those of high conductivity) 
to its mean diameter is the same for all the shells of the concentric as- 
sembly. In this case the shielding factor for two high conductivity 


shells is 
4 4 


where g equals the square of the ratio of the inside radius to the out- 
side radius for both high and low conductivity shells, 7 is equal to 
1—g and uz is the ratio of high to low conductivity of the alternating 
shells. 

A convenient formula for the shielding factor of multiple concen- 
tric shells has been developed by Walker and Sterne (6, 7). Assuming 
all the high permeable shells to be of the same thickness, ¢,, and of 
approximately the same mean radius, R, and all the low permeable 
shells to be of thickness, ¢,, two constants e; and é2 are defined as 
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Other constants are defined as 


= I — — + 


a 
B = é2, 
Y = 
6=1— — é2. 


The shielding factor is then 
F P(x:)"*) + Q(x2)"*}, 
where x; and 4 are the roots of the equation 
x? — (a + 6)x + (ad — By) = 0. 
P and Q are found from a solution of the two simultaneous equations 
Px, + Qxe = 1, 
Px? + Qx2? = 1/2a+B+7+ 5), 


and m is the number of high conductivity shells. 


Fic. 3. Simple thermostated oven. 


The discussion of shielding against temperature gradients has so 
far been concerned only with metal cylinders in the absence of sources 
of heat. In the construction of a gravity-meter thermal oven, it is de- 
sirable to wind the heating coil on one of the high conductivity cham- 
bers. Such a case is represented in Fig. 3 where a cylindrical metal 
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box, closed at both ends and having a wall d cm thick, is surrounded 
by an insulating material of thickness r. The ratio of thermal conduc- 
tivities is again equal to u. Assume that the thermostat is located in the 
wall of the cylinder a distance h/2 from either end. Also assume that 
a temperature gradient exists on the outside of the chamber parallel 
to the axis of the cylinder and equal to A¢/h where Ag is the difference 
in ambient temperature between the two ends. If a uniform winding 
of resistance wire covers the area of the metal cylinder, the thermostat 
will so regulate the flow of current in the wire that the heat generated 
per unit area will be constant everywhere on the surface of the cylinder 
and its end plates. If this constant is denoted by c’ watts/cm? the 
action of the thermostat provides that 
C= 


where TJ is the inside temperature and ¢ is the ambient temperature 
opposite the thermostat. The loss of heat per unit area as a function of 
x, the distance along the axis of the cylinder, is 


The net production of heat per unit area is then c’—b, which is 
simply (—A¢/r)(x/h—4). Because the wall of the cylinder is thin com- 
pared to its radius and length, it is possible to approximate the rate 
of heat production per unit volume as (—Ad@/rd) («/h—4) and applica- 
tion of the continuity equation equates this expression to —yV°T. 
For this one-dimensional case of heat flow parallel to the wall, 


a?T ( x 
dx? yurd\h 


dT Ad (x? 


where N is the constant of integration and is approximately equal to 
—1/u times the longitudinal flow of heat per unit area from the end 
plate of the cylinder (x=0). To a first approximation the end plate is 
at uniform temperature (rh) and the flow of heat from it to the 
cylinder wall can be taken as zr*(c’ —b,~0) = (Ad/27)xr? which corre- 
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sponds to a flow per unit area into the end of the cylinder of Agr/4rd. 
Substituting for N we have, 


dT Ad (= x “) 


dx wyrd\2h~ 2 4 
and 
rh 
urd \12 4 


where now, AT is the difference in temperature between the ends of 
the cylinder caused by a corresponding difference in ambient tempera- 
ture of Ad. 

Although the use of alternating layers of high and low conductive 
material is a logical way to obtain good temperature distribution, it 
leads to expensive construction and excessive weight. To economize 
on the weight of metal oven walls the writer has resorted to the use 
of several thermostats spaced symmetrically around a single oven. 
Each thermostat controls the production of heat in its own immediate 
neighborhood and even though each thermostat may operate at a 
slightly different temperature the resulting temperature distribution 
within the oven remains more nearly constant than is possible with a 
single thermostat. 

The best that can be expected of any thermostat is that it hold its 
own temperature constant and if other points are held constant it is 
only because they are in good thermal contact with the thermostat. 
With plural points of control there are no regions of the oven remotely 
located relative to a controlled point. 

Plural points of control seem to be superior even to the use of a 
thermostat, such as a resistance bridge, that can be spread uniformly 
over the oven wall. Although the resistance bridge averages the tem- 
perature of the oven and regulates the heat supply accordingly, it is 
obviously possible for the temperature distribution within the oven to 
change without changing the average temperature. 


ESSENTIAL COMPONENTS OF GRAVIMETER OVEN . 


The features so far presented serve as a practical guide in the de- 
sign of a thermal oven for a gravimeter. 

Since very close control of the temperature is necessary it is very 
probable that self sustained oscillations will exist. Highest sensitivities 
are obtained when the smallest temperature change which a given 
thermostat can detect, is used to switch the heat source entirely on or 
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off; i.e., use of the highest possible values of s jn continuous control 
degenerates into relay-type thermostats. It is therefore desirable to 
provide an attenuating layer (a layer of material for which pc/k is 
large) between the thermostat and the inside of the oven. See Fig. 4. 

To insure a uniform temperature within the oven a distributing 
layer should surround the attenuating layer. For highest efficiency 
this distributing layer should be made of aluminum containing a plu- 
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Fic. 4. Doubly thermostated oven wall. 


rality of thermostats, or it should be divided into two or more layers 
separated by insulating material. To economize on total heat consump- 
tion the heating coils and thermostat should be located on the inner- 
most distributing layer. The thermostat should be in as good thermal 
contact as possible with both the heat source and distributing layer 
in order that xo be small. The resulting high frequency for the thermo- 
stat cycle materially improves the efficiency of the attenuation layer. 

Mercury-in-glass thermostats have been found very satisfactory 
for use with gravimeter ovens because of their stability and the 
simplicity of the auxiliary control equipment which they require. 
When their sensitivity is so high that the backlash is of the order of 
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.02°C. they are sensitive to vibration; i.e., road shocks cause “‘chatter- 
ing” of the mercury column against the contact electrode. The chatter- 
ing is caused by the fact that the thin glass bulb containing the heavy 
mercury deforms when accelerated. This deformation can be greatly 
reduced by inserting the mercury bulb in a steel tube filled with mer- 
cury. The mercury bath cancels the inertia of the mercury in the bulb 
and also provides a good thermal connection between the bulb and its 
environment. Introduction of time delay in the control circuit is un- 
desirable because it increases the over- and undershooting in the ther- 
mostat. 

If the glass bulb of the thermostat is too thin the mean operating 
temperature will vary slightly with barometric pressure; a negligible 
effect except in areas of large variations in altitude. 

An insulation layer surrounds the distributing layer. The insula- 
tion should of course have a low value of k to economize on the heat 
supply as well as to increase the value of u and improve the efficiency 
of shielding against temperature gradients. However, in addition to a 
low k, it is important that the insulation have as large a value of 
pc/k as possible in order to attenuate short-time fluctations in am- 
bient temperature. Thus wool felt is preferable to kapok even though 
both have the same conductivity. 

To reduce long-time fluctuations in the ambient temperature it is 
convenient to surround the insulating layer with a second thermo- 
stated distributing layer. This second distributing layer is held at an 
approximately constant temperature several degrees lower than the 
inside temperature. Reference to Eq. 18 shows how the resulting 
constancy of (IT—@) and corresponding smallness of d@ will aid 
greatly in the reduction of dT to a very small value. Voltage regula- 
tion is also desirable to reduce dV /V toa low value. This second heating 
and distributing layer acting with its inner companion also gives the 
oven the advantage of multiple shielding against temperature gradi- 
ents. 

Finally, a second layer of insulation is added and it should have the 
same properties as the first insulating layer. On first consideration 
one may be inclined to make this second insulating layer quite thick 
in order to economize on the amount of heat necessary to maintain the 
outer heating layer at its proper temperature. However, the tempera- 
ture difference between the inner and outer heating layers requires 
that a certain amount of heat, q’, flow from the former to the latter 
and then continue outward through. the second insulating layer. The 


| | 


168 DAYTON H. CLEWELL 


heat q’ flowing through a thick outer layer of insulation would produce 
a large temperature drop and it would be necessary to limit the rise in 
-ambient temperature to a value rather low relative to the inside 
temperature. 

Even in the case of a single heating layer it is not desirable to use 
too much insulation because, as the ambient temperature approaches 
the inside temperature, the heating cycle becomes so unsymmetrical 
as well as of such low frequency that the attenuating layer is no longer 
effective. If the inside temperature is raised so as-to be well above the 
highest ambient temperature the corresponding increase in heat sup- 
ply that is then required, cancels the economy of thick insulation. 

A cross-section of what might be considered an ideal oven wall for 
a gravimeter is shown in Fig. 4. The heavy line is a graph of the 
mean temperature through the wall and the shaded area the tempera- 
ture variation. 

The writer wishes to thank Mr. H. C. Cortes and Dr. J. P. Minton 
of the Magnolia Petroleum Company for their kind permission to pre- 
sent this discussion. 
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THE ANALYTIC BASIS OF GRAVITY 
INTERPRETATION* 


D. S. HUGHESt{ 


ABSTRACT 


One method of gravity interpretation involves the use of analytic continuation 
processes. In this discussion the resolving power of this method is tested numerically. 

Using hypothetical structures comprising single and double blocks, a surface- 
gravity profile is derived. Using these values as an “observed gravity” profile, the 
“continuation” method is applied to compute the gravity at intermediate depths. Com- 
paring these computed values with the actual (directly computed) gravity profile at 
these depth-planes, the resolving power of the continuation method is demonstrated. It 
is shown that a very high precision in the observed data is necessary for very accurate 
resolution of structures. 


The problem of deriving the underlying structure from observa- 
tions of the gravitational attraction at the surface has received consid- 
erable attention. A large number of graphical and mechanical schemes 
of computation have been evolved for calculating the gravitational 
attraction of a given mass distribution, but the more difficult inverse 
problem of computing the mass distribution from the observations 
has received no satisfactory treatment. This latter problem is in fact 
known to have, in general, no unique solution. One common method 
of obtaining a solution is to compute the gravitational attraction of 
various simple mass configurations and then compare these computed 
values with the observed values. By varying the depth and size of the 
assumed configurations it will usually be possible to reproduce the 
observations with any desired accuracy. While this process is of great 
practical utility it must be clearly recognized that it does not give a 
unique solution and the actual structure that gives rise to the obser- 
vations may bear little or no resemblance to the structure derived by 
this method. 

Several writers have called attention to the fact that the solutions 
are more determinate for shallow than deep structure. Shallow and 
deep in this connection are terms relative to the size of the structure. 
In order to gain a practical advantage from this fact various attempts 
have been made to find a method of computing the gravitational at- 
traction at any given depth below the surface from observations made 
at the surface. Thus the object of this method is to take the surface 
gravity map and translate it into the map that would be obtained 
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were it possible to make the observations at some fixed depth below 
the surface. Obviously this depth may be chosen so that the map is ob- 
tained on a plane as near the structure as desired. Within certain 
limitations this problem may be solved and thus, theoretically, from 
a surface gravity map we should be able, by constructing maps on 
successively deeper planes, to attain any desired degree of resolution. 
Actually this method suffers from the same limitations as the method 
of assuming simple mass distributions and then varying them in such 
a way as to fit the observations. In fact a modification of this latter 
process can be shown to give a rigorous solution of the problem of com- 
puting gravity at depth from the observed values at the surface. 

In computing the anomalies from the observations, corrections are 
usually made for the shape of the earth, elevation, topography, isos- 
tasy, etc., and thus the principal features of the earth are accounted 
for and the computed anomaly is effectively the attraction of anoma- 
lous density distributions located in an outer shell of the earth 100 
kilometers or less in thickness. This shell is spheroidal but as its thick- 
ness is small compared with its radius a restricted neighborhood may 
be considered plane. Its upper surface (actual surface of the earth) 
will be taken as the (x, y)-plane and the z-axis will be taken positive 
upward. 

Then g(x, y, z) is the attraction of anomalous masses located below 
the plane z=o at every point (x, y, z). In particular g(x, y, o) is the ob- 
served anomaly at the point (x, y). H. M. Evjen' has shown that from 
g(x, y, ©) it is possible by a simple integration to compute g(x, y, z) for 
z greater than zero and furthermore that the successive derivatives of 
g with respect to z may be computed at all points for which z=0. With 
these derivatives it is then possible to compute g(x, y, 2) for values of 
z<o by means of Taylor’s expansion. This is essentially the standard 
process of analytic continuations. 

If we denote by g(x, y,—) the gravitation attraction on a plane h 
below the surface, the attraction at the surface is 


te ete g(u,v, — h)hdudv 


Since g(x, y, 0) is known and g(u, v, —h) unknown this expression is an 
integral equation for the unknown g(u, v, —h). In the following we 
shall be concerned only with one set of values of (x, y) and hence this set 


1H. M. Evjen: “Place of Vertical Gradient’? Gropnysics, January, 1936. 
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may without loss of generality be taken to be (0, 0) and only the z 
coordinate will be specified explicitly. 

Polar coordinates r, 0, are more convenient than the rectangular 
coordinates u, v. With these alterations Eq. (1) becomes, 


g(r, 0, — h)hrdr dé 


In Eq. (2) the variable @ is involved only in g(r, 0, —h). Hence, setting 


g(r, — = f(r, 0, — mao, 


Eq. (2) may be written 


g(r, — h)hrd 


This equation may be solved by the method of iteration or repeated 
approximation. As a first approximation we may assume that 2(r, —h) 
in the integral is identical with Z(r, 0) and then compute an approxi- 
mate value of gravity at the surface. 

Denoting the degree of approximation by subscripts, 


&0(— h) = g(o) 

g(r, o)hrdr ; 
g0(0) = (2+ (4) 

The computed gravity go(o) at the surface is in error by 

— h) — 2(r, 0) |h rdr 
(r? h?)3!2 
This error may be regarded as a new anomaly to be accounted for and 
as before we may assume that a better approximation will be given by 


g(o) — go(o) = (5) 


g(o) — go(o) = gi(— — g(0), 
or solving, 


gi(— h) = 2g(0) — go(o). (6) 


But it is obvious by reference to Eq. (4) that 


80(0) g(h), 
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or that go(o) is the anomaly that would have been observed on a plane 
h above the surface and 


gi(— A) = 2g(0) — g(h). (7) 


The substitution of this improved approximation for g(r, —/) in Eq. 
(3) will give a better approximation for the surface gravity. 


— h)hrdr 
-f (r2 + h?)3/2 (8) 


Substituting gi(—/) from Eq. (7) this equation may be written 


= 2g(h) — g(2h), (9) 
and the error of this approximation is 
g(o) — gi(o) = g(o) — 2g(h) + g(2h) (10) 


0 
and again we may assume the approximate value of the error to be 
— gi(o) = g2(— h) — 
or substituting from equations (7) and (9) 
82(— h) = 3g(0) — 3g(h) + g(2h). (11) 


This approximation may in turn be substituted in Eq. (3) and 
g2(o) computed. 


= 3g(h) — 3g(2h) + g(3h). (12) 
By continually separating this process it is found that the mth approxi- 
mation to g(—h) is 
n+1 
k=0 k+1 
and the mth approximation to the observed anomalies is 
n+1 
k=0 k+1 


Thus far it is not known whether or not this process is convergent, 
that is whether or not the successive approximations are really better 
approximations, In the form given above it is not possible to answer 
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this question. The convergence of the process depends upon the g’s 
and cannot be investigated until they are specified. Eq. (14), it will 
be noted, is obtained from (13) merely by adding / to the argument 
of all g’s. Hence it is sufficient to investigate only Eq. (13). 

From Eq. (3) it is obvious that 


g(r, o)kh rdr 


By a well known relation? 
kh 
+ 


0 
So that we may write 
g(kh) o)rar f o(Ar)AdX, 
0 0 


or 


g(kh) = o)Jo(Ar)rdr. (17) 


The first form is not convergent for k=o, although obviously it is for 
all positive values of k. The second form is convergent for h=o, and 
is merely the Fourier-Bessel integral. Substituting Eq. (17) for the g’s 
in Eq. (13) there is obtained, 


; “alr, o)Jo(Ar)rdr; (18) 
but 
and 


gn(— h) = foe [1 — (1 — o)Jo(Ar)rdr. (19) 


2A. G. Webster: “Partial Differential Equations of Mathematical Physics.” 
Stechert & Co., 1927, page 406. 
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In this equation if m becomes large without limit 


g(— h) = f f “ir, o)Jo(Ar)rdr. (20a) 


If g(x, y,o) is constant in the y direction or the observed values are 
not a function of y, Eq. (13) still applies; however, Eq. (202) becomes, 


g(— h) = fet, cos Ax dx. (20b) 


The equations are rigorous formulae for computing the attrac- 
tion on a plane at depth / from observations at the surface. The in- 
tegrals will be convergent only if / is less than the depth to the shal- 
lowest singularity of g(r, 0). Eqs. (20) are of little or no practical 
utility in the actual solution of the interpretation problem. 

If g(r, o) were given in analytical form it would offer a solution of 
the problem but the observations are point to point numerical values 
and the numerical solutions of Eqs. (20) are hopeless. These equations 
do, however, offer confirmation of the limiting process and of the finite 
forms given by Eggs. (13) and (14). These latter equations must oper- 
ate under the same limitations as Eqs. (20), i.e., 4 must be less than 
the depth of the shallowest singularity of g(r, o). A priori, this depth 
is not known; however, if for some selected # the successive approxi- 
mations from Eq. (13) do not approach a limit, or if the integrals in 
Eqs. (20) do not converge, the selected / is too large and a smaller 
value must be chosen. . 

In contrast to Eqs. (20), numerical computations may be carried 
out with Eq. (13). It is only necessary to compute g(kh) for a selected 
number of planes above the surface. Thus at any point g,(—/) may 
be computed from Eq. (13). This process, while long, is not impossible. 
Eq. (13) is, however, the condition that the (w+1)th difference of a 
function given at the (n+ 2) points —h, 0, h, 2h, . . . nh, shall be zero; 
in other words, that the (7+2) values can be fitted by a polynomial 
of the nth degree. It is thus simply the polynomial extrapolation for- 
mula. 

It is well known that large errors are inherent in such an extrapola- 
tion and are unavoidable. These errors are not introduced specifically 
by the above mathematical derivation, but are inherent in any scheme 
of interpretation. 

The mathematical derivation of Eqs. (13) and (20) is merely an 
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analytic formulation of the method of interpretation by assuming a 
primary distribution and then introducing corrections so that the final 
structure will reproduce the observations. It may be urged that in the 
present case the first approximation was very crude and that a better 
first approximation would give a better convergence. This, however, is 
not true. A good first approximation makes fewer.steps necessary to 
achieve a selected accuracy but does not increase the accuracy of the 
final result. This is easily seen to be true since if the process converges 
at all some one of the successive g,(—/) will approach as closely as we 
please to the true value of g(—z). This particular g,(—h) may be 
regarded as the first solution and the process continued from there, 
the final accuracy will not be changed. 

A very useful result which may be derived from these equations is 
the effect of errors in observations. Error in this connection denotes 
any deviation of the observed value from the true attraction of the 
subsurface structure. Thus not only instrumental errors are included 
but also errors due to improper or incorrect reduction formulae and 
errors due to assumptions of incorrect densities for the strata above 
sea level. Ordinarily the topographic correction is not computed and 
this introduces an additional error. 

If Eq. (13) is used to compute the anomaly g,(—/) below the sur- 
face g(o), the observed value directly above this point is multiplied 
by n, the degree of approximation, and any error in g(o) is thus in- 
creased m times. The g(#) of Eq. (13) is multiplied by m(m—1)/2 and 
consequently its error is increased even more; on the other hand, g(/) 
is computed by an integration over all the surface observations and 
thus random errors will, to a large extent, be cancelled. The degree of 
cancellation cannot be formulated exactly but it is apparent that if the 
station spacing is small compared with 4, the cancellation of random 
errors will be much more effective than for spacings of the same magni- 
tude as h. If the error in all the g(kh) was the same as that in g(o), the 
probable error of gn,(—/) would be the probable error of g(o) multi- 
plied by the square root of the sum of the squares of binomial co- 
efficients in Eq. (13); on the other hand, if the errors of g(kh) were 
zero, the error in g,(—/) would be times the error of g(o). Thus for 
a fifth order approximation, the error will be multiplied by some factor 
greater than 5 and less than 16. 

Hence the degree of approximation that can be used is determined 
by the accuracy of the data and the allowable error in final results. 
Eqs. (20a) and (20b) correspond to an infinite degree of approximation 
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and hence can only be applied to data of infinite accuracy, or in other 
words, to data given by an analytic formula. 

As an example of the application of these equations the two struc- 
tures shown in Fig. 1 have been chosen. The computed gravity profiles 
for each structure for several] different depths are shown. A density 
difference of 0.3 has been assumed in computing the profiles. For the 
structures actually drawn a constant should be added to the curves for 
the broader structure since its base level is different. However, in the 
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Fic. 1. Computed gravity profiles for two hypothetical structures at several as- 
sumed depths, Z. Dashed lines are for single block structure; solid lines are for double- 
block structure. 


reduction of field observations we never obtain an absolute anomaly 
but only a relative one. In other words, the actual numerical anomaly 
has no meaning. The interpretation is based on the excess or defi- 
ciency of the anomaly at one point as compared with another. These 
two structures shown here have been chosen in such a way as to be 
different as possible and yet to give as similar anomalies as possible. 
It will be noted that if the double structure reaches the surface and the 
single structure is 500 feet deep the two g curves are quite different. 
However, if both structures are placed 1000 feet deeper the g errors are 
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very similar showing a maximum difference of only 0.3 milligals. At 
a depth of 2000 feet the difference in the g curves is only 0.15 milligals. 
and at 3000 feet it is less than 0.10 milligals. 

In order to demonstrate the use of the equations derived earlier, 
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Fic. 2. Showing the resolution obtainable by “‘continuation methods” applied to 
gravity profiles. Circular and triangular points respectively, show “continuation” of 
observed gravity (Z= 2000) for double and single block structure of Fig. 1 at Z= 1000 
and Z=o planes. Compare points with actual gravity in these planes as shown by the 
curves. 


the computed g curves at 2000 feet above the structure have been 
selected. These curves and also the curves for 1000 feet and zero are 
shown in Fig. 2 plotted to a larger vertical scale. The upper curves 
(z= 2000 feet) have been used as the observed curves and by means of 
Eqs. (13) and (15) the observations at 1000 feet and o have been com- 
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puted. These points are shown by circles for the double structure and 


triangles for the single structure. 

Fairly satisfactory reproductions of the curves at 1000 feet are 
obtained and the resolution is increased somewhat. The computed 
points at z=o show a resolution but so far off the actual curves that 
it is evident that the corresponding structures would be quite different 
from the original structures. 

A fifth degree approximation was used in computing these points. 
A more accurate reproduction cannot be obtained by taking a higher 
order formula. The computation of the original curves is accurate to 
0.01 mgl. To increase the accuracy of the project points it would be 
necessary to use not only a higher order approximation but also a more 
accurate original curve. The departure of the projected points from 
the curves should not be regarded as an error. These points represent 
curves at depth that will reproduce the original curves within the 
limits of accuracy of the observations. In other words, these projected 
points represent possible interpretations of the observations within 
the accuracy of the observations. 
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THE CROSBYTON ANOMALY, SOUTHEASTERN 
CROSBY COUNTY, TEXAS* 


ETHEL WARD McLEMORE,}+ PAUL WEAVER,{ anp 
DONALD C. BARTON$ 


FOREWORD 


The purpose of this study has been the interpretation of the most 
probable structure producing the slightly subcircular anomaly at 
Crosbyton, as mapped by the pendulum, torsion balance, magneto- 
meter, and dip needle. Work was begun in 1935 under the direction of 
Dr. Barton, who was to have written the geological considerations. 
All calculations were completed in 1938, but unfortunately the geologi- 
cal interpretation was not written by Dr. Barton before his death, 
July 8, 1939. At my request, since he had discussed the Crosbyton 
anomaly with Dr. Barton many times, and was quite familiar with the 
work we were doing, Mr. Weaver kindly consented to write the geo- 
logical interpretation and help me with the completion of the paper. 
It is due to his knowledge and able assistance that the paper has been 
brought up to date in the light of subsequent drilling by the Gulf Oil 
Corporation. 

HISTORY AND LOCATION 


The Crosbyton high lies nearly in the center of the high plains of 
West Texas, approximately 45 miles southeast of Lubbock, and near 
the common corner of Crosby, Dickens, Kent, and Garza Counties. 
The high is named for the County seat, Crosbyton, of Crosby County, 
which is about fifteen miles north from the apex of the geophysical 
picture. 

The average elevation of the area is approximately 2300 feet, with 
low relief. Drainage is to the southeast through the White River, one 
of the tributaries of the Brazos. There are no particular difficulties in 
transportation, nor in surface conditions which obstruct geophysical 
operations. 

Although until recently there has been no important oil or gas 
development near the Crosbyton high, it marked a milestone in the 
economic geology of Texas, because it was from a well drilled at Spur, 


* Presented at the Eleventh Annual Meeting, Houston, Texas, April 2, 1941. 
} Research Geophysicist, Dallas, Texas. 

t Chief Geophysicist, Gulf Oil Corporation, Houston, Texas. 

§ Died July 8, 1939, Houston, Texas. 
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in Dickens County, on the northeast side of the feature, that the first 
potash salts were reported from the Permian. Dr. J. A. Udden reported 
on this occurrence,! and it led to the drilling of a number of special 
test wells, subsequent to the World War which demonstrated the im- 
portance of a domestic source of potash in this country. This campaign 
of test wells was instrumental in the development of the present oper- 
ating potash mines in New Mexico. 

The Spur well was a local effort, drilled under pioneer conditions 
(Swenson fee #1, NW/4, Sec. 225, H. & G. N. RR., Dickens County), 
and fifty-two months elapsed from spudding until the total depth of 
4489 feet was reached in November, 1913. It was at that time the 
deepest well that had been drilled in Texas. 

In January, 1926, drilling was begun on a well in Garza County 
by the Deep Rock Company, about four and three-fourths miles south 
and two and one-fourth miles west from the northeast corner of Garza 
County. It was not until October 23, of that same year that the total 
depth of 3502 feet was reached (Deep Rock, Swenson #1, Center Sec. 
36, Blk. 2). 

In 1923 the first large production in West Texas was brought in at 
Big Lake, about 100 miles south of Crosbyton, causing immediate 
attention by geologists. This was almost simultaneous with the intro- 
duction of geophysical instruments in the Gulf Coast; but, since these 
were presumed to be useful primarily for the discovery of salt domes, 
they were not taken to West Texas immediately. However, since 
granite was known to be present in West Texas in the Panhandle, and 
since large areas seemed to be difficult to work as to surface geology, 
it was natural that attention should be given to an attempt to find 
some relation of buried granite features to structure in areas where 
surface geology was inconclusive. Thus the dip needle was tried in 
1925, in the Panhandle. In 1926 the Askania magnetometer came into 
use in a general reconnaissance of the entire Permian Basin by several 
companies. It was on the basis of detailed work with this instrument 
that the Gulf Production Company drilled the #1A Swenson. There 
also at this time the refraction seismograph and the torsion balance 
were introduced. ; 

The first mention in the literature of the Crosbyton high was in 
February, 1928, at the meeting of the American Institute of Mining 
and Metallurgical Engineers in New York City, where Dr. Barton 


1 J. A. Udden, Deep Boring at Spur, University of Texas Bulletin, 363, October 5, 
1914. 
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presented a paper.” In the discussion Mr. DeGolyer® asked about the 
Crosbyton high; and the following comments by the senior author are 
inserted her, since the original publication is not commonly available.* 


It seems to me it must be a big batholith. Mathematically, it could 
be accounted for by a large structural high of limestone, rising to 
within 2500 feet of the surface and sloping off rather gently to a depth 
of three miles, at a distance of about twenty miles from the crest... . 
We can get a series of pictures of possible structures that are all 
equally probable mathematically. In general they may have a mini- 
mum depth or a maximum depth. In this case we could not bring the 
structure within 2500 feet of the surface by any reasonable assump- 
tion, but the top of it might be 5000 feet just as easily. The Gulf 
Production Company has drilled a test below 5000 feet without finding 
anything but the normal stratigraphic section. 


The well by the Gulf Production Company to which the senior 
author referred was located in the northeast corner of Garza County 
(NW corner, Sec. 67, H. & G. N. RR., Blk. 2) and was the Swenson 
#1A, completed March 24, 1928, at 5105 feet. On zonal correlation 
this well was correlated as thirty-six feet lower than that of Deep Rock. 

Subsequently, in 1938, another well (#1B Swenson) has been drilled 
by the Gulf Oil Corporation, which found granite at 8085 feet, con- 
tinuing in that formation to the total depth of 8104 feet; and followed 
by a second well (#1C Swenson) to a depth of about 8414. All the cal- 
culations in this article were completed prior to the drilling of these 
wells. In the geologic interpretation which concludes this paper will 
be found comments on comparison of the calculations with the results 
of these wells, the only deep tests drilled in the general area since the 
geophysical work was done. 

The following is an abstract of the projected paper written by Dr. 
Barton in 1937: 


A large slightly subcircular anomaly has been mapped by pendu- 
lum, torsion balance, magnetometer, and dip needle, near Crosbyton 
in Northwest Central Texas. The width of the anomaly is approxi- 
mately 4o miles, its gravitational amplitude 45 milligals and its 
magnetic amplitude 2500 gammas. A well which was drilled to a 
depth of 5105 feet, near the crest of the anomaly, is carried as normal 
by the subsurface geologists. According to mathematical calculations, 
the anomaly must be the effect of a flattened stock like mass whose 
vertical dimension is several miles. From geological considerations it 
must be the effect of a plutonic rock. 


2 D.C. Barton: “The E*tvi's Torsion Balance Method of Mapping Geologic Struc- 
ture,”’ Geophysical Prospecting, 1929, A.I.M.E., pp. 416-466. 
3 Op. cit. p. 467. 
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GEOPHYSICAL INTERPRETATION 
Data from independent surveys by four different instruments were 
used in an attempt to work out a composite profile across the anomaly, 
from which was calculated the probable structure producing this 
anomaly. These surveys were: pendulum determinations made by the 


Fic. 1. Compromise isogams between torsion balance and Bouguer pendulum 
anomalies. Insert—Observed torsion balance gradients. 


United States Coast & Geodetic Survey during the fall of 1933; torsion 
balance survey made by D. C. Barton for the Amerada Petroleum 
Corporation in May, 1927; magnetometer surveys made by Amerada 
with the Gepege Vertical Variometer early in 1927, and by the Texas 
Company with the Askania Vertical Field Balance about the same 
time; dip needle survey made by Amerada with the Hotchkiss Super 


Dip Needle in 1927. 
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Helmert’s formula of 1901, g = 978.030(1+0.005302 sin? 0.000007 
sin? 20) was used in the computation of theoretical values of gravity at 
the twelve stations, where gravity observations and latitude deter- 
minations were made by the United States Coast & Geodetic Survey. 
The observed values were adjusted to a common latitude of 33°30’ 
and sea-level; applying the elevation correction according to the 
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Fic. 2. Magnetic isogams. Insert—Deflection of dip needle. 


Bouguer formula. From these calculated theoretical values of gravity 
and the adjusted observed values were derived the isostatic anomalies 
(Zobs.--Ztheor.) Which are shown on the map in Fig. 1. 

Maps of the four surveys were plotted on a common scale of one 
inch equals to 10,000 feet (Figs. 1 & 2). From each of these a corre- 
sponding east-west profile was calculated. These are shown on Fig. 3. 
From these profiles a compromise profile was derived (Fig. 3a). Since 
the details of the slope of this profile are principally from torsion bal- 
ance data, and its azimuth is nearly east-west, it is designated as the 
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“observed gradient profile” of the Crosbyton anomaly, or a U,, profile 
in the language of the torsion balance. 

The graphic method of calculation of gradient profiles from dif- 
ferent assumed structures was carried through, employing a U,, chart 
of 9:10:10 proportions,‘ for twenty stations selected at equal distances 
along the east-west cross-section (Figs. 3 and 3a show location of 
stations from which the values of the “‘observed gradient profile” could 


Fic. 3a. Observed gravity-gradient profile. 


be read for comparisons with each assumed structure). In order to 
have the complete structure come within the dimensions of the chart, 
the scale of one inch equal to 2500 feet was used, the “observed gra- 
dient profile’ being replotted on this scale. All results were subse- 
quently replotted on the scale of one inch equal to 10,000 feet. 

The U,, chart was used for computation of three-dimensional 
figures because it is obvious from the geophysical contours that the 
Crosbyton high is due to some mass of nearly equal dimensions in any 


‘D. C. Barton: “Calculations in the Interpretation of Observations by the Eotvos 
Torsion Balance,’”’ Geophysical Prospecting, 1929. Amer. Inst. Min. Met. Eng., 1929, 
Pp. 480-504. 
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horizontal direction. In the cross-sections comparing computed and 
observed gradients, only two dimensions are shown, except in Fig. 10 
where the best fit between calculated and observed values indicated 
that the mass has the figure of a parallelopiped inclined to the horizon- 
tal as to its major dimension. In Fig. 10 the cross-section is shown in 
perspective. 

The first assumed structure was a finite rectangular block 
(A’B’C’'D’) with axis median with respect to the maximum and mini- 


wer 


SYMMETRICAL RECTANGLE 
DENSITY +020 


d 


Fic. 4. Symmetrical rectangle—with profiles. The gravity profile is calculated 
for a three-dimensional mass with square horizontal cross-section. 


mum of the anomaly. With variations of depth and position the gra- 
dient values at twenty stations were calculated, with constant size of 
the block. To obtain the average density contrast of the assumed struc- 
ture, the average cf the ratio of the observed to the calculated values at 
each station was taken; and, since the density enters in the chart as 
a multiplier, this average gives the density contrast. The least square 
method was used in selecting the most probable depth and position of 
this assumed structure. The two calculated gradient profiles of the 
two positions of the structure A’B’C’D’ which gave the least sums of 
the squares of the deviations were compared with the “observed 
gradient profile” (Figs. 4 & 5). A comparison of these two profiles with 
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the observed at each station shows that in the center of the cross- 
section there is a fair coincidence, with a much faster increase of the 
observed values eastward, and a slow tendency toward constancy of 
the two calculated profiles. Westward there is a consistency as far as 
Station — VII, when the observed values increase rapidly, while the 
calculated values seem to approach constancy, by slow decrease. 

As this seemed to point to the possibility of a composite structure 
of a finite block of one density and a block of another density and same 


OENSITY +0 20 


Fic. 5. Symmetrical rectangle; second position—with profiles. The gravity profile 
is calculated for a three-dimensional mass with square horizontal cross-section. 


cross-section extending downward from the first to infinity, the values 
of the infinite block MNOP, with top at the base of the other block, 
were calculated. This composite structure and the gradient profile 
produced are shown in Fig. 6, with a comparison with the “observed 
gradient profile.”” While the inclusion of the infinite block gives a 
better fit eastward, positive values much higher than the observed 
throw it altogether out of comparison around Stations O to —VII 
inclusive, though giving the same increase of values past Station— VII 
as that of the observed. 
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Fic. 6. Finite block with square horizontal cross-section 
together with block MNOP 
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Fics. 7 & 8. Structure “X”’ with profiles. 
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In an attempt to get a better idea of the incongruency of the as- 
sumed structures, calculations for stations farther out from the ob- 
served gradient profile were located, both east and west. Calculations 
of the values of gravity at these stations proved that the tendency to 
decrease slowly westward, approaching zero as a limit, and increase 
slowly eastward continued with the outer stations. 

Another structure ‘““X”’ was assumed, a parallelopiped instead of 
a rectangle, dipping 20° westward, and deeper than the best position 


SECOND MOST PROBABLE STRUCTURE 
DENSITY* 020 


Fic. 9. Second most probable structure—with profiles. Block QRST has square 
horizontal cross-section. 


calculated for the rectangle (Fig. 7). Comparison of the profile of this 
structure with the observed profile shows it a little too far eastward, 
and dipping at too large an angle (Fig. 8). 

Two other parallelopiped structures were assumed: QRST dipping 
at an approximate angle of 223° westward, and W’X’Y’Z’ dipping at 
an approximate angle of 15° westward, and moved over considerably 
westward of the center of the anomaly (Figs. 9 and 10). The sums of 
the squares of the deviations from the observed values of these calcula- 
tions were the smallest yet obtained. Their profiles are shown com- 
pared to the observed in Figs. 9 and 10. Except for an unaccountable 
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discrepancy at Station +II, QRST, with a density contrast of 0.2, 
seemed to be the most probable structure. However, considering this 
discrepancy, W’X'Y’Z’, with an average density contrast of 0.23, gives 
a profile of nearest coincidence with the observed (Fig. 10). 

Neither a composite profile of W’X’Y’Z’ with the infinite block, 
nor variations of the size and shape of the infinite block seemed to 
give a better profile to the observed. 


Fic. 10. Most probable structure; in perspective—with profiles. 


Although these assumed masses as to shape, size, depth, and posi- 
tion are the only ones discussed here, many variations were calculated. 
Since many of the readers will be interested in the outline of the calcu- 
lations of one of these profiles, and the numerical comparison of ob- 
served and computed gradient values, the summary computation table 
for the mass finally selected as the best fit (Fig. 10) is included. 


GEOLOGICAL INTERPRETATION 


The combination of a magnetic high and a gravity high in an area 
where the sediments were already known to have density at least the 
average of sedimentary rocks, suggests a basement rock of plutonic 
source, and that it was either a large monadnock or an intrusion into 
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SUMMARY COMPUTATION TABLE 
W 
Obs./Cal. 
Cal. d 
Density X0.23 (dev.) 
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Least Square Sum= 288.52 


the sedimentaries. The #1B Swenson found granite at a depth of 
— 5810 feet subsea, which was within 40 percent of the depth to the 
plutonic rock calculated on Fig. 10, at the point where the well was 
drilled. The #1C Swenson did not reach granite at 5936 feet 
subsea, showing that the mass is dipping to the west, as was calcu- 
lated on Fig. 10. 

It would not have been expected that granite could give such 
density contrast with the limestones and shales above it as is actually 
found, and this would suggest that the granite is not the main con- 
stituent of the plutonic mass, but that the mass is principally a more 
basic rock than granite. Since the #1B Swenson drilled only 19 feet in- 
to the mass, we have no confirmation of this deduction. ; 
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A THEORY OF THE DISTRIBUTION OF 
RADIOACTIVITY IN MARINE 
SEDIMENTARY ROCKS* 


PAUL WEAVERt 


ABSTRACT 


Gamma-ray measurements of marine sedimentary rocks indicate that in general 
the slower the rate of deposition the greater the radioactivity. The most recent meas- 
urements of the radioactivity of sea water suggest an increase with depth, which indi- 
cates a very slow settling of small particles. It is believed that a geologic time-scale of 
sedimentation can be established, subject to certain corrections due to fixation of radio- 
active minerals by organisms. 


Rays are emitted when there is a transformation from one radio- 
active element to another. The elements which possess this property 
naturally are potassium, rubidium, and the elements of the three 
radioactive series—uranium, actinium, and thorium. The rays are of 
three classes, alpha-, beta-, and gamma-, and for each element the 
rays are specific as to frequency of emission and as to penetrating 
power through different media, therefore the number of gamma-rays, 
for example, which arrive at a point inside a well is dependent upon 
the kind of element which emits the rays, the distance of the emittor- 
atom from the observation point, and the absorbing power of the 
medium between, which may be drill mud or other fluid in the well, 
casing, or geological formation. 

Certain gamma-rays are very penetrating, and therefore the kind 
of fluid in the well, presence or absence of casing, etc., reduce the pro- 
portion of the rays which arrive inside the hole only slightly, so that 
variation of their frequency from one point to another is proportional 
to the change in percentage, within the strata alongside the wall of the 
well, of a radioactive element or elements. 

Apparatus has been inserted in a number of wells within recent 
months, whereby, with a connection to a recording device at a surface, 
a record or “‘log’’ can be made of the variation with depth of the total 
of the gamma-rays which originate in the drill mud or other fluids in 
the well and of those which are emitted from the strata forming the 
wall of the well. The apparatus which has been used is either a (1) 
Geiger-Mueller Counter or (2) Ionization Chamber. The Geiger- 
Mueller Counter set-up was described in GEopuysics™. 

* Presented at the Annual Meeting, Houston, Texas, April 2, 1941. 


+ Gulf Oil Corporation, Houston, Texas. 
t Figures refer to references at end of article. 
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There are other methods which are being used to measure the pro- 
portion of radioactive elements in rock or soil, which are not readily 
adaptable to use for logging wells, and are principally applied in the 
laboratory to hand specimens. These determine one of the following: 


(1) Amount of alpha-rays from a gaseous radioactive element, 
either thoron or radon, or both from which the parent element 
is computed on equilibrium assumption. 

(2) Amount of helium present, from which is computed total 
emission during past geologic time on assumption that helium , 
is all derived from radioactive disintegration. 

(3) Ratio of various isotopes of lead, since some of these are end 
products of radioactive disintegration. 

(4) Velocity of beta-rays which are specific for each radioactive 
element. 

(5) Secondary radiations emitted by exposure to a radioactive 

source. This method has recently been applied to well-logging. 


Since there is such a spectrum of rays due to the many radioactive 
elements in each of the three families of uranium, actinium and — 
thorium, and since they are usually all present in many minerals found 
in sedimentary rocks, the values determined by the different methods 
are not convertible to those measured by the Geiger-Mueller Counter 
or the ionization chamber run into a well, nor are gamma-ray well logs 
precisely expressible in any unit to correct to a chemical formula ex- 
pressing the kinds and proportions of each radioactive element present 
in the strata. It will be the purpose of this paper, however, to show 
that practically the readings of these logs have a definite interpreta- 
tion, and to propose a theory that the variation in these readings from 
one kind of sedimentary rock has an application which gives the geolo- 
gist a new time-scale. 

It is already recognized that gamma-ray logs are useful to show 
the boundaries between different kinds of marine sedimentary rocks 
and for correlation between wells some distance apart.* The results of 
many wells are concordant in their pictures of the four great types of 
marine sediments. These are—(1) muds, clays and shales; (2) sands, 
sandstones and quartzites; (3) calcareous oozes, limestones, dolomites 
and chalks; (4) evaporites—salt, anhydrite and some limestones and 
dolomites. The relative gamma-ray intensity of these are as follows:— 


* In some publications, a log by the use of the ionization chambers has been called 
a “Radioactivity log.’? Gamma-ray is a more precise term, and also has priority (see 
Reference (1)). It is hoped that the commercial operators will return to the “Gamma- 
ray”’ to identify logs made by ionization chambers. 
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Evaporites lowest 

Sands low 

Clays high 

Limestones usually low, occasionally quite high. 


Where measurements are made from surface samples for correla- 
tion with wells, there will be discrepancies unless fresh samples are 
used. Surface work on radioactivity content where the soil processes 
have begun will generally give higher values than characteristic of the 
formation because weathering results in an increase in the more in- 
soluble rock minerals, among which are those with the radioactive 
elements as a constituent. 

The marine sedimentary deposits of the present oceans have been 
studied by bottom samples taken by dredging and by coring. The 
muds accumulate on the bottom by flocculation of kaolin and similar 
colloids by the sea water; the sands by mechanical transport due to 
wave action; the lime oozes by decrease in carbon dioxide in the water 
and by action of organisms; and the evaporites are produced by pre- 
cipitation when the ocean water has been concentrated by evaporation 
—anhydrite when 70% of the original water is evaporated, salt when 
90% has been evaporated. 

Geologists have speculated concerning the rate of sedimentation 
in the ocean, using various methods based on the rate of transport into 
the ocean by rivers and other physical factors. Some recent values are 
as follows, expressed in thickness after compaction :— 


Shales 1 foot per goo years® 
Sandstones 1 foot per 450 years 
Limestones 1 foot per 2250 years® 
Evaporites 1 foot per 3 years. 


There is no unanimity concerning the numerical values given 
above or even the order of magnitude, as the methods used by different 
geologists to calculate them have been based on different processes. 
However, there is agreement as to the general relation expressed 
above, that evaporites form very much faster than any other type of 
sediment; sandstones are faster in rate of deposition than shales which 
again are faster in rate of deposition than limestones.* 

* Limestones are the present form not only of what were originally calcareous 
oozes, but of evaporites, of clastics made of lime fragments, of coquina beds, and of 
secretions of organisms. Only as to the limestones formed from calcareous oozes, can 
we use the rate of deposition mentioned. Where other conditions are present, limestone 
source material may be deposited very rapidly. 
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It will be observed that the geologist has estimated a comparative 
time-scale for the deposition of the three classes of marine sediments— 
muds, sands and evaporites, according to which they are arranged in 
the same order that these classes respond to gamma-ray logging, that 
is, the faster the rate of deposition estimated by the geologist, the less 
the gamma-ray reading. Interestingly enough, this is outstanding in the 
case of halite, or common salt, the principal evaporite, which is 
the least radioactive of all rocks known—in fact, the best place in the 
world to calibrate a Geiger-Mueller counter is in a salt mine, provided 
the salt is pure or nearly so, as has been recognized before there was 
any gamma-ray logging. 


= Evaporites (Salt, etc) 
— Sandstones 

Shales 


Limestones (Poor) 
| (Rich) 


Arbitrary Scale (Linear) 


Fic. 1. Relative Gamma-Ray Activity of Sedimentary Rocks. 


There is much evidence in literature of measurements of the radio- 
activity of sedimentary rocks taken on hand specimens in the labora- 
tory and on soils overlying the rocks, taken in the field, but a recapitu- 
lation seems unnecessary, first, because the technique was not always 
standard and precise; and second, because the logging of thousands of 
feet of drilled wells is equivalent to many thousand times all the other 
data. However, the author is satisfied that the relative radioactivity 
of the four classes of sediments from all the evidence is substantially 
as shown in Fig.1.* 


* Sediments contain locally special mineral zones. For example, a sandstone de- 
posited near an old land mass of pre-Cambrian rocks will have much heavy mineralsand 
therefore may be exceptionally high. Phosphate nodules in shales and limestones may 
be radioactive due to replacements of base radicles by radium, etc. Such special condi- 
ditions are quite useful for geologic correlations, giving sharp, narrow peaks on the 
gamma-ray logs, but are not characteristic of thick zones. 
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Assuming an inverse correlation therefore between rapidity of 
sedimentation in the ocean and degree of radioactivity in the sedi- 
ments, there is naturally a query as to the cause. If we can find this 
cause in some physical process in the ocean water, we may be able to 
measure it in operation today and obtain a numerical conversion 
factor from gamma-ray reading to years per foot of deposition, after 
gamma-ray variations have been calibrated so that logs taken under 
different operating conditions are comparable. 

Now how do the radioactive elements occur in the ocean? Recent . 
measurements® show that the radioactive content of sea water is quite 
low, but increases with increasing depth at greater than a linear rate. 
It is also evident that in the deep quiet bottoms, far from land and 
where wave action is minimized, the oozes on the ocean bottom have 
very high radioactivity compared either to the average of all oceanic 
sediments or to the sea water in contact with them. The author pro- 
poses that these observations support a theory that the mechanism of 
accumulation of radioactive minerals is the sedimentation of very 
small discrete particles of the parent elements in mineral combination 
falling from the air through the sea water to the bottom at a rate which 
is statistically constant. A certain amount of these particles is dis- 
solved enroute, a larger amount is fixed temporarily by organisms 
which are strongly selective, such as kelp, but the total effect of all of 
these interruptions to the downward settling of the particles is con- 
sidered of second order, because of the very small solubility and of the 
relative paucity of the selective organisms. 

What is the source of the mineral particles assumed to be slowly 
settling? They can come only from four sources: (1) river muds of 
extreme fineness and stable in sea water; (2) volcanic dust from erup- 
tions; (3) meteor clouds; (4) general dust circulation in larger wind 
movements. River muds and volcanic dusts would be locally concen- 
trated and should result in areas of very high radioactive sediments 
irrespective of whether they be shales or sandstones or evaporites. 
Perhaps more widespread application of the gamma-ray log will show 
such areas, but they are not now known. Only general wind circulation 
and meteor clouds are widely enough distributed in time and space to 
give uniform conditions over large areas. Meteors were for many years 
considered non-radioactive, but recently it has been reported that 
some are. However, the total tonnage which falls annually into the 
atmosphere and reaches the earth’s surface unconsumed, is small, and 
the size of individual fragments probably too large for slow sedimenta- 
tion in the ocean. 
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Therefore there seems left as the source for the particles of radio- 
active material only the general air circulation carrying the finest par- 
ticles derived by wind action (aeria] erosion) on the land surface. This 
material would be of small diameter. It would contain principally 
those minerals of low solubility in water vapor, of considerable hard- 
ness, and of a greater proportion of heavy minerals than the average 
of the earth’s crust, because it was derived principally from residual] 
soils rather than from fresh rock. Such properties are those which 
would supply in relative abundance the zircon, monazite, etc., in 
which heavy elements are most frequent. 

If this dust from general air currents be the source of the radio- 
active minerals, its rate of supply to the surface of the ocean would be 
expected to be statistically constant over a geological unit of time. Its 
rate of sedimentation through the sea water would also be constant 
because biological data indicates constancy of composition, viscosity, 
temperature of sea water at different depths through long periods of 
time. 

The best check on this hypothesis would be the collection of large 
samples of atmospheric dust and their radioactive analysis, also a 
statistical investigation of grain size and mineral composition, so that 
a computation could be made of the rate of sedimentation and loss by 
solution in sea water. These figures compared with radon now in solu- 
tion in sea water would check the general theory for the uranium 
family. Similar measurements for the thorium family would be some- 
what more difficult, but not impossible. As far as present literature is’ 
concerned, the evidence does not seem to be forthcoming for a first 
numerical calculation. 

This theory is worth investigation for several reasons. In the first 
_ place, it would indicate that limestones are poor for correlation of 
geologic time, and that there are different limestones with different 
histories. The limestone reservoirs are important as oil producers, and 
are puzzling in their rapid areal variations (reefs). 

Again, when a shale thickens from one well to another, there is 
now uncertainty as to whether the difference is due to variation in 
amount of supply of the original mud, or whether there has been down- 
warping so that more deposition time was involved where the shale is 
thicker. The establishment of downwarping as the cause of the thick- 
ening would help the search for stratigraphic traps, if used in connec- 
tion with isopach maps made by the reflection seismograph. 

Finally, the establishment of a geologic time-scale for sediments will 
help in studies of compaction and regional metamorphism, especially 
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in monoclines and geosynclines, which are the tectonic provinces of 
most of our oil fields. There might even be cases where the paleonto- 
logical correlations could be checked over some distance. 
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DISCUSSIONS AND COMMUNICATIONS 


WAR EFFORTS ADVISORY COMMITTEE 


At the annual meeting of the Society in Denver on April 24, 1942, a resolution 
proposed by B. B. Weatherby was unanimously adopted by those present as follows: 

Resolved, That the Society of Exploration Geophysicists cooperate fully in the 
War Effort of our country by aiding in every possible way the maintenance of an ade- 
quate exploration program and also further the placement of all available geophysical 
personnel in suitable positions with the armed forces and other governmental agencies; 
also that the Society render assistance to all geophysical laboratories and shops desiring 
to do work for the government and for these purposes a committee be appointed by 
the President. 

To put into effect this resolution the President, Dr. Peacock, immediately ap- 
pointed a committee to be known as the War Efforts Advisory Committee. Those 
appointed were the following: B. B. Weatherby—Chairman, Henry C. Cortes, R. D. 
Wyckoff, Herbert Hoover, Jr., W. G. Green and F. Goldstone. 

At the time of going to press, it has not been possible to ascertain that all of these 
appointees will be able to serve but it is anticipated that a committee substantially as 
above will be functioning shortly. Accordingly, members who have problems which 
they believe fall within the scope of the activities of this committee should communicate 
with the chairman. 


COMMITTEE ON EXPLORATION UNDER O.P.C. 


On March 9g, under the Production Committee for District 3, as defined by the 
Office of the Petroleum Coordinator, a subcommittee was organized, to be known as 
the Committee on Exploration, comprised of Messrs. Sidney A. Judson, Phil F. Martyn, 
R. R. Morse, Ed. W. Owen, Paul Weaver, and F. H. Lahee, chairman. This Committee 
on Exploration is a fact-gathering group of men, not representing any particular com- 
pany association, or special interests, but serving as representatives of the petroleum 
industry. The function of the Committee is to assemble, as soon as possible, data bear- 
ing (1) on the progress of wildcatting (i.e., exploration for oil by drilling outside of known 
producing pools, or below the deepest known pay sands in producing pools); (2) on 
prospects for petroleum not yet exploited; (3) on methods of exploration to be used in 
search for new prospects; (4) on the quantity of steel needed in exploring for new pros- 
pects; and (5) on the quantity of steel needed for wildcat drilling for oil on new pros- 
pects. 

An early and comprehensive study of these and related data is vital in the program 
of coordinating the petroleum industry for national defense. Therefore, it is expected 
that all persons and organizations approached for information toward this end will be 
ready to cooperate. The committee will welcome any constructive suggestions. 


F. H. LaAwEE, Chairman 
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ELECTRICAL PROSPECTING 


U.S. No. 2,268,106, L. W. Blau, Iss. 12/30/41, Appl. 5/13/39. Assign. Standard Oil 

Development Co. 

Radio Wave Prospecting—A method of locating subsurface structures, particularly 
faults, by recording and mapping the field strength or combinations of field strengths 
of one or more distant radio transmitters. 

U.S. No. 2,270,325, N. W. Matthews, Iss. 1/20/42, Appl. 1/18/41. 

Auxiliary Electrode for Ground Resistance Measurement—A ground electrode having 
at its upper end two collars between which a weighted sliding handle may be used 
either to drive the electrode into the ground or to drive it out after use. 


GEOCHEMICAL PROSPECTING 
U.S. No. 2,257,170, L. G. Howell, Iss. 9/30/41, Appl. 12/27/34-1/21/38. Assign. 
Standard Oil Development Co. 
Method for Locating Hydrocarbon Deposits in the Earth—A method for determining 
hydrocarbon gases in the soil by sucking soil gas from a shallow bore hole, passing the 
gas through a liquid which dissolves hydrocarbons, and then analyzing the liquid. 


U.S. No. 2,269,889, L. W. Blau, Iss. 1/13/42, Appl. 2/27/39. Assign. Standard Oil 

Development Co. 

Process for Locating Valuable Subterranean Deposits—A method of prospecting in 
which soil samples are tested with a reagent which produces a color change when ap- 
plied to soil containing hydrocarbon consuming bacteria which have been subjected to 
action of hydrocarbons. 

U.S. No. 2,270,299, L. Horvitz, Iss. 1/20/42, Appl. 3/23/39. Assign. E. E. Rosaire. 

Geochemical Prospecting—A method of soil analysis involving treatment of soil 
samples with hydrochloric or other acid to dissolve carbonates and release trapped 
hydrocarbons. 

U.S. No. 2,272,645, E. E. Rosaire, Iss. 2/10/42, Appl. 11/21/41. 


Geochemical Prospecting—A method of geochemical prospecting in which shallow 
holes are treated with acid and the released gas pumped off and analyzed. 


GRAVIMETRIC PROSPECTING 


U.S. No. 2,262,165, D. H. Clewell, Iss. 11/11/41, Appl. 7/3/40. Assign. Socony- 
Vacuum Oil Co. 


Gravity Meter—A null reading gravity meter in which the main spring is displaced 
transversely by an adjustable secondary spring so as to change the lever arm of the 
main spring. 

U.S. No. 2,264,342, D. Silverman and J. L. Bible, Iss. 12/2/41, Appl. 5/12/39. Assign. 

Stanolind Oil and Gas Co. 


Geophysical Instrument Mounting—A gravimeter mounting which eliminates re- 
moval of instrument from truck and affords protection in transit. 


* Abstracts by O. F. Ritzmann, Gulf Research & Development Co. 
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SEISMOGRAPH PROSPECTING 


U.S. No. 2,257,187, J. E. Owen, Iss. 9/30/41, Appl. 12/21/39. Assign. Geophysical 

Research Corporation. 

Seismic Surveying—A method of eliminating ground roll from a seismic record 
which consists in combining with the output of the usual velocity type vertical geophone 
either the output of a displacement type horizontal instrument or an acceleration type 
horizontal instrument. 


U.S. No. 2,259,478, C. G. Morgan, Iss. 10/21/41, Appl. 10/1/38. Assign. Consolidated 
Engineering Corporation. 
Geophysical Exploration System—A method of seismic prospecting in which con- 
tinuous control is obtained by shooting into a spread from regularly spaced shot points 
more distant than those adjacent to the spread. 


U.S. No. 2,261,321, P. S. Williams, Iss. 11/4/41, Appl. 6/17/38. Assign. Standard Oil 
Development Co. 
Seismic Exploration Method—A method to facilitate the identification of reflections 
on a seismic record by recording along with the customary traces one or more traces 
showing average rectified amplitude. 


U.S. No. 2,263,519, O. F. Ritamann, Iss. 11/18/41, Appl. 12/2/38. Assign. Gulf 
Research & Development Co. 
Seismograph Prospecting A pparatus—Eliminating seismic interference falling within 
the useful range of frequencies by superimposing on the usual band pass filter a rela- 
tively sharper band absorption filter. 


U.S. No. 2,264,596, E. M. Shook and R. W. Olson, Iss. 12/2/41, Appl. 10/23/40. 
Assign. Socony-Vacuum Oil Co. 
Geophysical Instrument—Method of making seismic shot moment or other time 
breaks identifiable through static by having shooting device start a relaxation oscillator 
thus repeating the break at regular known intervals. 


U.S. No. 2,265,513, K. E. Burg, Iss. 12/9/41, Appl. 6/10/39. Assign. Geophysical 
Service, Inc. 
Prospecting Method and A pparatus—A method of seismic prospecting in which each 
detector is equipped with a radio transmitter and recording is accomplished via radio 
instead of cables. 


U.S. No. 2,265,538, J. P. Minton, Iss. 12/9/41, Appl. 7/9/38. Assign. Socony-Vacuum 
Oil Co. 
Gain Control for Seismograph Amplifiers—A manual gain adjusting circuit having 
a linear relation between gain and position of control index. 


U.S. No. 2,265,768, L. F. Athy, H. R. Prescott and R. F. Hughes, Iss. 12/9/41, Appl. 

7/20/40. Assign. Continental Oil Co. 

Method of Logging Bore Holes—A seismic method of well logging making use of a 
continuous wave transmitter to set up standing waves and measuring the impedance 
changes as the transmitter is lowered down the well, also exploring the standing wave 
field with several vibration receptors lowered into the well with the transmitter. 
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U.S. No. 2,266,040, H. Hoover, Jr., Iss. 12/16/41, Appl. 9/16/37. Assign. Consolidated 
Engineering Corp. 


Geophysical Prospecting Receptor Circuits—A seismic detector-amplifier input cir- 
cuit having a shielded isolation transformer whose primary is balanced to ground and 
whose secondary may be interconnected to other channels to permit mixing. 


U.S. No. 2,266,041, H. Hoover, Jr., Iss. 12/16/41, Appl. 4/18/38. Assign. Consolidated 

Engineering Corp. 

Geophysical Prospecting Receptor Circuits—A seismic detector-amplifier input cir- 
cuit having two shielded transformers in cascade, the primary being balanced to ground 
and the intermediate circuits being interconnected and adjustable to eliminate cross- 
feed and permit controlled mixing. 


U.S. No. 2,266,837, H. F. Wiley, Iss. 12/23/41, Appl. 7/24/39. Assign, Consolidated 
Engineering Corp. 
Method and Apparatus for Translating Seismic Waves—A geophone whose weight 
and base area are so proportioned that the frequency of maximum response when 
planted is above the frequency range of reflections. 


U.S. No. 2,267,357, O. F. Ritzmann, Iss. 12/23/41, Appl. 1/16/40. Assign. Gulf 
Research & Development Co. 
Apparatus for and Method of Seismograph Recording—A method of mixing seismo- 
graph channels optically by having the recorded light beam reflected successively from 
two galvanometers. 


U.S. No. 2,267,858, C. H. Dix, Iss. 12/30/41, Appl. 9/8/41. Assign. Socony-Vacuum 
Oil Co. 
Method of Seismic Surveying—A method of refraction profiling which requires 
shooting in only one direction. 


U.S. No. 2,268,130, M. M. Slotnick, Iss. 12/30/41, Appl. 2/8/40. Assign. Standard 
Oil Development Co. 
Method of Geophysical Investigation—A method of outlining a subsurface structure 
by placing a seismic detector in a well penetrating the structure and shooting into it 
from shot points on the surface. 


U.S. No. 2,269,453, E. L. Gayhart, U.S.N., Iss. 1/13/42, Appl. 1/12/38. 


Device for Detecting Dis placements—A displacement meter or seismometer in which 
a coil moves in an alternating magnetic field so that the alternating current induced is 
proportional to the displacement. 


U.S. No. 2,269,890, L. W. Blau, Iss. 1/13/42, Appl. 9/30/39. Assign. Standard Oil 
Development Co, 
Seismic-Electric Pros pecting—Seismic-electric detectors consisting of two current 
electrodes and one or more pairs of potential electrodes at stations where it is desired 
to pick up seismic disturbances. 


U.S. No. 2,271,864, P. M. Honnell and L. W. Dickerson, Iss. 2/3/42, Appl. 6/10/38. 
Assign. The Texas Co. 
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Seismic Wave Detector—A suspended coil electromagnetic geophone provided with 
electromagnetic and electrostatic shielding. 


U.S. No. 2,272,741, T. E. Failing, Iss. 2/10/42, Appl. 2/28/38. Assign. George E. 
Failing Supply Co. 
Apparatus for Seismic Surveying—A non-metallic non-resonant shot hole casing 
which eliminates transmission of shock up the casing. 


U.S. No. 2,272,984, O. F. Ritzmann, Iss. 2/10/42, Appl. 7/21/39. Assign. Gulf Re- 
search & Development Co. 
Seismograph—A geophone having a band pass response characteristic which con- 
sists of two suspended masses with different resonant frequencies and a transducer 
whose output depends on the difference of their individual response characteristics. 


WELL LOGGING AND SURVEYING 


U.S. No. 2,261,563, F. Rieber, Iss. 11/4/41, Appl. 2/18/39. 

System for Measuring Earth Conductivity—A method of electric well logging while 
drilling by passing an alternating current between the drill stem and a distant surface 
electrode and measuring by a null method the magnetic field created at a third point 
on the surface. 


U.S. No. 2,262,419, L. F. Athy and H. R. Prescott, Iss. 11/11/41, Appl. 4/9/38. 
Assign. Continental Oil Co. 
Method of Electromagnetic Logging—A well logging method using surface current 
electrodes on each side of the well and a nondirectional horizontal magnetic field de- 
tector lowered into the well. 


U.S. No. 2,265,978, D. E. Batchelder, Iss. 12/16/41, Appl. 10/4/39. Assign. Lane- 
Wells Co. 
Alternating Current Electro-Logging of Well Bores—An A. C. well logging method in 
which simultaneous curves are drawn of the in-phase potential, quadrature potential 
and phase angle. 


U.S. No. 2,266,071, R. G. Piety, Iss. 12/16/41, Appl. 5/31/39. Assign. Phillips Petro- 
leum Co. 
Well Surveying Device—A method of electric well logging using a surface electrode 
and a very long wire electrode in the well and measuring the current distribution leaving 
the wire by interposing a sliding shielding electrode connected to a meter. 


U.S. No. 2,268,137, H. M. Evjen, Iss. 12/30/41, Appl. 12/11/39. Assign. Shell Devel- 
opment Co. 
Electrical Well Logging System—A four electrode method of logging uncased holes 
which makes a sequence of A.C. and I).C. measurements controlled by synchronous 
commutators so as to obtain three log curves simultaneously. 


U.S. No. 2,268,138, H. M. Evjen, Iss. 12/30/41, Appl. 12/11/39. Assign. Shell Develop- 
ment Co. 
Electrical Well Logging System—A four electrode method of logging uncased holes 
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which simuitaneously gives the limited and extended radius resistivity (A.C.) curve and 
also the spontaneous (D.C.) potential curve. 


U.S. No. 2,268,256, W. S. Knouse, Iss. 12/30/41, Appl. 2/3/40. Assign. 1/2 to National 

Lead Co. 

Apparatus for Surveying Deep Wells—A geometrical well surveying instrument 
using an opaque liquid to change the amount of light falling on three photo cells for 
determining inclination and a shaded compass card to change the amount of light 
falling on a photo cell for determining direction. 


U.S. No. 2,268,627, D. Silverman, Iss. 1/6/42, Appl. 10/6/39. Assign. Stanolind Oil 
and Gas Co. 

Well Logging—A method of well logging using a drilling fluid under high pressure, 
substituting a liquid having a different physical or chemical property, releasing the 
pressure to allow drilling fluid to reenter the well, and measuring the chemical or 
physical property at various levels in the well. 


U.S. No. 2,268,682, E. R. Webb, Iss. 1/6/42, Appl. 10/9/37. Assign. Sperry-Sun Well 
Surveying Co. 
Well Surveying Instrument—An inclinometer using a pendulum and recording by 
means of an electric spark. 


U.S. No. 2,260,269, F. W. Jessen, Iss. 1/6/42, Appl. 7/15/40. Assign. Standard Oil 
Development Co. 
Well Logging—A method of locating breaks or ends in casing which consists in 
lowering a galvanic electrode into the well and measuring the potential between this 
electrode and the casing. 


U.S. No. 2,260,573, A. Barry, Iss. 1/13/42, Appl. 4/20/30. Assign. Standard Oil 
Development Co. 
Adjustable Measuring Wheel—A well logging cable or wire line measuring wheel 
having adjustable circumference. 


U.S. No. 2,260,717, D. Hering, Iss. 1/13/42, Appl. 12/13/40. Assign. Sperry-Sun Well 
Surveying Co. 
Well Surveying Device—An inclinometer in which colored ink marks a cylindrical 
record so arranged that the ink need not be removed between runs. 


U.S. No. 2,271,051, J. M. Pearson and G. A. Smith, Iss. 2/3/42, Appl. 2/12/36, 

Assign. Sperry-Sun Well Surveying Co. 

Electrical Prospecting Method and Apparatus—A well logging apparatus which 
records in the hole the effect of surrounding formations on a high frequency oscillating 
circuit. 

U.S. No. 2,273,215, J. Neufeld, Iss. 2/17/42, Appl. 4/22/40. Assign. Well Surveys, Inc. 

Well Surveving Method and Apparatus—A method of determining the point of 


entry of formation fluid into a cased hole by making repeated radioactivity logs under 
reduced hydrostatic pressure. 


PATENTS 205 


U.S. No. 2,273,363, E. Lipson, Iss. 2/17/42, Appl. 11/28/39. 

Method for Electrical Investigation of Cased Drill Holes—A four electrode method 
of electrical well logging using two electrodes on the surface and using current values 
which are disproportionate to the potential, e.g., o.2 to 1.0 amp. 


MISCELLANEOUS 


U.S. No. 2,268,808, J. O. Davis, Iss. 1/6/42, Appl. 12/6/39. Assign. Henry Hughes & 
Son, Ltd. 
Recording Device—A stylus type recorder for echo sounding having devices for 
changing speed, zero correction and phasing adjustments. 


U.S. No. 2,260,414, J. P. Woods, Iss. 1/6/42, Appl. 10/16/39. Assign. Shell Develop- 
ment Co. 
Multiple Recording Galzonometer—A permanent magnet oscillograph unit having 
elements mounted in separately removable cartridges and using a condensing lens to 
concentrate light on the galvanometer mirrors. 
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CONTRIBUTORS 


Dayton H. CLEWELL received his B.S. 
degree in 1933 and his Ph.D. degree in 1936, 
both from the Massachusetts Institute of 
Technology. He was employed by the C. K. 
Williams Company of Easton, Penna., from 
1935 to 1938 as a physicist working on optical 
problems associated with pigments. Since 1938 
he has been associated with the Geophysical 
Laboratory of the Magnolia Petroleum Com- 
pany where he has been concerned mainly with 
the development of gravimeters. He is the 
author of three papers in the field of experi- 
mental physics. He is a member of the Ameri- 
can Physical Society and of Sigma Xi. 


Dayton H, CLEWELL 


D. S. Hucues received his A.B. degree 
from the University of Kentucky in 1926 and 
his Ph.D. degree from the University of Chi- 
cago in 1931. From 1931 to 1933 he held a 
National Research Council Fellowship at the 
California Institute of Technology. He has been 
engaged in geophysical work since 1935 when 
he became associated with the Shell Oil Com- 
pany, by whom he is still employed. Dr. 
Hughes’ work has been primarily in the field of 
gravity interpretation. 

He is a member of the American Physical 
Society and the Society of Exploration Geo- 
physicists. 


D. S. HuGHEs 
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ETHEL Warp McLemore (Mrs. R. H. 
McLemore) received her B.A. degree from 
Mississippi Woman’s College in 1928, majors 
in chemistry and mathematics; M.A. degree 
from University of North Carolina in 1929 in 
chemistry and mathematics; work on Ph.D. 
at University of Chicago in 1931. She entered 
the Research Department of Geophysics of 
Humble Oil & Refining Company in 1934, 
having taught chemistry and mathematics at 
Mississippi State College for Women in 1930- 
32. She remained with the Humble until 1937. 
Since that time she has worked on mathe- 
matical interpretation of geophysical data, par- 
ticularly as pertaining to gravimetrical instru- 
ments and the correlations of their results, and 
private teaching of advanced mathematics and 
French. At present Mrs. McLemore is working Era. Wan Mclaucen 
on a collection of unpublished papers by Dr. 

D. C. Barton, to be published posthumously by the American Association of Petroleum 
Geologists. She is a member of Sigma Xi Society and the American Association of 
University Women. 


Raymonp G. Prety received his B.S. de- 
gree from the University of Chicago in 1932 
and his Ph.D. in Physics from Yale in 1936. He 
was engaged in sound motion picture develop- 
ment for Warner Brothers Pictures, Inc. during 
1930 and 1931 and again during 1935-and 1936. 
He joined the research staff of the Phillips 
Petroleum Company in 1937 where he is still 
employed. His work as section chief in charge 
of geophysical research has been largely con- 
cerned with problems of instrument develop- 
ment. He is a member of the American Physical 
Society and the Society of Exploration Geo- 
physicists. 


RayMOND G. PIEtTy 
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A. SHARPE 


CONTRIBUTORS 


JosEpuH A. SHARPE received his B.S. degree 
in physics from the University of Arizona in 
1928, and his Ph.D. degree in physics from the 
University of Wisconsin in 1934, after which 
he spent one year at the Massachusetts Insti- 
tute of Technology as a General Education 
Board (Rockefeller Foundation) fellow in geo- 
physics. He has been associated with the joint 
geophysical laboratory of the Stanolind Oil and 
Gas Company and the Western Geophysical 
Company since 1935, first as seismologist, and, 
since 1937, as chief physicist. Dr. Sharpe’s work 
has been primarily concerned with the initia- 
tion, propagation, and characteristics of elastic 
waves, although he has also worked on seismic 
and geomagnetic instrumentation and inter- 
pretation, and on gravity interpretation. He is 
a member of the American Physical Society, 
the American Mathematical Society, the Amer- 
s, the Society of Exploration Geophysicists, the 


American Geophysical Union, the Seismological Society of America, and the Tulsa 


Geological Society. 


NEIL R. Sparks received his B.S. degree 
in Electrical Engineering from the University 
of California, Berkeley, in 1932. He continued 
his graduate studies in the field of Seismology 
until 1933 at which time he entered the employ 
of the U. S. Coast and Geodetic Survey in their 
program of earthquake research. Since 1936 he 
has been associated with the Western Geo- 
physical Company and the Stanolind Oil and 
Gas Company in their combined Geophysical 
Laboratories. Mr. Sparks’ work has been pri- 
marily the interpretation of seismic data. He 
is a member of the Seismological Society of 
America and the Society of Exploration Geo- 
physicists. 


NEIL R. SPARKS 
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REVIEW 


Geophysics in War, by C. A. Heiland, Professor of Geophysics, Colorado School of 
Mines. Published in Vol. 37, No. 1 of the Colorado School of Mines Quarterly. 
Price $1.00 postpaid. 


The science of geophysics is concerned with the physical phenomena in and the 
transmission of energy through the earth, the water, and the air. Its applications in 
war and preparation for war fall into two groups: military operations, and the location 
of essential minerals. 

In the combat zone acoustic triangulation or sound ranging helps to locate hostile 
guns. Listening devices determine the approach of submarines or airplanes. Sound- 
ranging methods are also employed in establishing the range of friendly artillery. Land 
and sea mines may be controlled by radio and acoustic impulses. Buried munition 
dumps, shells, and bombs can be located by the radio detection devices. Communication 
between battlefleet vessels is often maintained by sonic and supersonic transmission. 
Position finding at sea is facilitated by radio-acoustic ranging. In periods of poor visi- 
bility vessels may navigate by echo depth sounding. A similar application is terrain- 
clearance determination by high-frequency radio from a plane. 

Planning of fortifications (gun emplacements, shelters, harbors) and the location 
of construction materials for them will be aided by seismic refraction, electrical re- 
sistivity, and dynamic ground testing. The same methods are used in all kinds of trans- 
portation problems involving the construction of railroads, highways, bridges, tunnels, 
etc., and in the planning of munitions plants. For the last, added protection may be 
afforded by static-ground-resistance investigations. 

Meteorology, an important part of geophysics, is applied in weather forecasting for 
military operations, maneuvers, in sound ranging, and acoustic airplane spotting. Data 
on rainfall, humidity, and flood frequency are useful in the planning of certain defense 
plants. Salvage operations, particularly the location of shipwrecks and practice weapons, 
are aided by echo-sounding and radio location methods. 

The second group of geophysical applications is concerned with the location of 
mineral resources, comprising, water, fuels, and strategic minerals. Water may be 
located under favorable conditions by electrical and seismic geophysical methods; salt- 
water contamination may be traced, and water wells may be tested by electrical logging. 
Geophysical methods mentioned before in connection with foundation investigations 
are applicable to water-control problems arising in connection with irrigation, flood- 
control, and power projects. Of vital importance in modern warfare is an adequate 
supply of petroleum. Geophysical exploration by magnetic, gravimetric, seismic-re- 
flection, and electrical well-logging methods occupies a prominent place in the oil in- 
dustry. Coals and lignite deposits may be mapped by geophysical methods. Magnetic, 
electrical, gravimetric, and seismic exploration methods are now widely used in a 
systematic, government-sponsored exploration program to uncover vitally needed de- 
posits of bauxite, chromite, manganese, mercury, nickel, tin, and tungsten. 
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MEMBERSHIP APPLICATIONS APPROVED FOR PUBLICATION 


The Executive Committee has approved for publication the names of the following 
candidates for membership in the Society. This publication does not constitute an 
election, but places the names before the membership at large. If any member has 
information bearing on the qualifications of these nominees, he should send it to the 
Secretary within thirty days. (Names of sponsors are placed beneath the name of each 
nominee.) 

ACTIVE 


Robert Frank Aldredge 

G. H. Westby, T. A. Manhart, A. J. Barthelmes 
Gerald Allan Burton 

H. R. Thornburgh, W. Hafner, Harry L. Thomsen 
Cecil Edmund Chesher 

W. Hafner, H. R. Thornburgh, Harry L. Thomsen 
C. Calhoun Cody, Jr. 

G. H. Westby, T. A. Manhart, A. J. Barthelmes 
Jesse Robert Cornett 

G. H. Westby, T. A. Manhart, A. J. Barthelmes 
James Edward Dorris 

G. H. Westby, T. A. Manhart, A. J. Barthelmes 
Laurence Meade Hubby 

George D. Mitchell, Jr., C. C. Zimmerman, Roy L. Lay 
Phillip John Kempsky 

W. T. Born, L. Y. Faust, B. B. Weatherby 
Theodore Koulomzine 

(Article ITI-D-1 of the Constitution) 
W. Warren Longley 

(Article III-D-1 of the Constitution) 
John Joseph Lynch S. J. 

James B. Macelwane, S. J., Louis B. Slichter, Cornelius G. Dahm 
John William Marston 

G. H. Westby, T. A. Manhart, A. J. Barthelmes 
Haydon Wilburn McDonnold 

George D. Mitchell, Jr., Roy L. Lay, C. C. Zimmerman 
Jacob Neufeld 

W. G. Green, T. A. Manhart, Alexander Wolf 
Albert Clark Reid 

Joseph A. Sharpe, Daniel Silverman, Jack R. Cooper 
George Hendricks Speirs 

G. H. Westby, T. A. Manhart, A. J. Barthelmes 
Edward Banks Stover, Jr. 

Joseph A. Sharpe, Jack R. Cooper, Neil R. Sparks 
Joel Howard Swartz 

J. Brian Eby, James B. Macelwane, S. J., C. E. Van Orstrand 
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John Gordon Cole 

Cecil H. Green, Phil P. Gaby, Frederick Romberg 
Larkin Earl Deacon 

F. Goldstone, Willy Hafner, M. O. Gibson 
Stuart Sawyer Giezendanner 

George D. Mitchell, Jr., Roy L. Lay, C. C. Zimmerman 
Otto George Holekamp 

F. Goldstone, W. Hafner, M. O. Gibson 
Ford Lawrence Johnson 

H. W. Rose, W. G. Smiley, Jr., L. G. Ellis 
Henry Leonhard Rase, Jr. 

Paul F. Clement, Herbert J. Ferber, Victor J. Pezdek 
Rudolph Gust Sohlberg 

F. E. Romberg, C. H. Green, K. E. Burg 
David McIntosh Steel 

T. R. Shugart, E. E. Wilson, H. Hallock Brown 
Thomas Scott West 

Cyril B. Smith, F. Goldstone, David S. Muzzey, Jr. 


TRANSFER TO ACTIVE 


W. T. Lea 
J. B. Meitzen 
Edward C. Reagor 
Harry H. Sisson 
ASSOCIATE 
Robert Wiley Adams 
Joseph LaConte, John W. Mathews, T. P. Ellsworth 


STUDENT 


Ralph Theodore Duffner 

W. S. Levings, R. Maurice Tripp, Dart Wantland 
Benjamin Franklin Howell, Jr. 

Beno Gutenberg 
Kennedy Robert Marble 

Dart Wantland 


Louis Charles Pakiser, Jr. 
Dart Wantland 
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PERSONAL ITEMS 


Dwicut E. Warp, formerly with the Carter Oil Company, has accepted foreign 
employment with the International Petroleum Company, P. O. Box 803, Guayaquil, 
Ecuador, S. A. 


GeEoRGE H. WESSLER, JR., associate member of the Society with N. K. P. M. in 
the Netherlands East Indies upon the outbreak of the Pacific war, is reported to have 
joined the U. S. Army. 


C. T. RANKIN should be addressed c/o Geophysical Research Corporation, P. O. 
Box 2040, Tulsa, Oklahoma. 


RoBERT BaRLING, formerly with the Western Geophysical Company, is now Chief 
Geologist and Engineer for the Miller Geophysical Company of Los Hagen, California. 
He is located at 830 18th St., Bakersfield, California. 


D. G. Ecan has returned from the Netherlands East Indies and is located in Crow- 
ley, Louisiana, with the mailing address P. O. Box 366. 


Horace S. Moore is serving as Field Draftsman and Instrumentman with the 
Southwestern Engineering Company of Los Angeles. He may be addressed at 2204 
Floral Drive, Whittier, California. 


Epwarp S. Foster, Jr., Well Logging Operator with the Halliburton Oil Well 
Cementing Company, is located at 604 Madison Street, Lafayette, Louisiana. 


C. H. Drx is now Acting Chief Seismologist with the United Geophysical Company, 
1255 East Green Street, Pasadena, California. 


Matcoim S. Morrison, Geophysicist with the Gulf Research and Development 
Company, should be addressed at Box 1342, R. D. 1, Frankstown Road, Wilkinsburg, 
Pennsylvania. 


KENNETH E. Bure has resigned from Geophysical Service, Inc. to accept a position 
as Technical Supervisor with the Stanolind Oil and Gas Company, Philcade Building, 
Tulsa, Oklahoma. 


GeEorGE A. GarrETT, formerly with the Independent Exploration Company, is 
acting as Research Associate in the Radiation Laboratory of the Massachusetts Insti- 
tute of Technology in Cambridge. His mailing address is 81 Elm Street, Belmont, 
Massachusetts. 


STANLEY W. Wi1cox, of Seismograph Service Corporation, should be addressed at 
709 Kennedy Building, Tulsa, Oklahoma. 


WALTER D. Barrp, Computer with Geophysical Service, Inc., is now located at 
P. O. Box 391, Fresno, California. 


J. M. KiLgorn, Magnetometer Party Manager with the Gulf Research and Devel- 
opment Company, is located in Lubbock, Texas, with the mailing address P. O. Box 
1436. 
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News Letter No. 4 of the U. S. Waterways Experiment Station, P. O. Box 631, 
Vicksburg, Mississippi, lists ““Subject-Headings for Material in Engineer Department 
Research Centers Library.”’ It is noted that printed matter, slides, and films may be 
obtained on a thirty-day loan basis by individuals connected with reputable firms or 
institutions. 


E. W. Jounson is acting as Seismograph Observer for the Cities Service Oil Com- 
pany, Bartlesville, Oklahoma. 


L. B. SLICHTER is serving on the National Defense Research Committee of the Office 
of Scientific Research and Development, roth Floor, 172 Fulton Street, New York, 


Joun W. TExTER has returned from foreign employment and is doing seismograph 
interpretation with the Gulf Research and Development Company, P. O. Drawer 
2038, Pittsburgh, Pennsylvania. 


Lr. G. W. Hostetter, formerly with the Petty Geophysical Engineering Company, 
is located in the Signal Corps General Developmental Laboratories, Fort Monmouth, 


Red Bank, New Jersey. 


AustIN N. STANTON is President of Petroleum Science, Inc., with offices in the Kirby 
Building, Dallas, Texas. 


Vincent J. Mercter, of the Lane-Wells Company, has the mailing address P. O. 
Box 221, Seminole, Oklahoma. 


W. E. Wa ttts has accepted foreign employment with the West India Oil Company, 
Apartado XXIII, San Jose, Costa Rica. 


On March 14th, the Rice Chapter of Sigma Xi invited Houston physicists and 
engineers to attend the annual lecture, presented by Dr. Percy W. Bridgman, Hollis 
Professor of Mathematics and Natural Philosophy of Harvard University. The subject 
was “Some Recent Work in the Field of High Pressures.” 


GeorGE B. Somers has advised the Society of his transfer to the Naval Bureau of 
Ordnance. As he is continually on the move, correspondence shou!d be addressed to his 
home, 4336 Bowser Avenue, Dallas, Texas. 


GeEorGE N. MEADE is serving as Jr. Seismograph Observer for the Magnolia Petro- 
leum Company, Box goo, Dallas, Texas. 


R. W. LrnpsEy, Geophysicist with the Gulf Research and Development Company, 
is located at P. O. Box 213, New Orleans, Louisiana. 


Lr. R. G. Paterson, R. C. N. V. R., should be addressed c/o Naval Officer in 
Charge, Vancouver, B. C., Canada. 


H. Wayne Hoytman, of the Gulf Research and Development Company, is on for- 
eign duty at Puerto Cabezas, Nicaragua. 


H. M. Zenor has left the Humble Oil and Refining Company to accept employment 
with the Sperry Gyroscope Company, Inc., Clinton Road and Stewart Avenue, Garden 


City, New York. 


; 
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P. M. KonkEL, of the Ohio Oil Company, has been transferred from Rochester, 
Indiana, to Marshall, Illinois. 


RAYMOND LAWRENCE SARGENT resigned from the Amerada Petroleum Corporation 
on March 23rd to accept a position as Chief Computer on Crew 20 of the National Geo- 
physical Company, Tower Petroleum Building, Dallas, Texas. 


CorporRAL CHARLES C. Bates, formerly with the Carter Oil Company, is on duty 
with Battery A, 2nd F. A. Obsn. Bn., Camp Bowie, Texas, as Chief Computer in the 
Sound Plotting Squad. 


DonutL HI1tIs is now serving as Geologist with the Capital Company, 1 Powell 
Street, San Francisco, California. 


F. F. CampBELL, Party Chief with the Geophysical Research Corporation, is lo- 
cated at 2416 Arbor Street, Houston, Texas. 


Joun C. Horrman, Party Chief with the United Geophysical Company, should be 
addressed at 1255 East Green Street, Pasadena, California. 


F. E. Curnutt returned from South American employment several months ago, 
and is now located at P. O. Box 11, Kentwood, Louisiana. 


Bulletin No. 14 of the U. S. Waterways Experiment Station Soil Mechanics Bulle- 
tin contains an abstract of the English translation of ‘Relations Between the Behavior 
of Clay Minerals with Respect to Water and the Physical Properties of Cohesive Soils,”’ 
by K. Endell, W. Loos, H. Meischeider, and V. Berg. A copy is on hand at Society 
headquarters. 


Francis A. Roserts, of the Carter Oil Company, is located at P. O. Box 349, 
Hazlehurst, Mississippi. 


Offices of the Mott-Smith Corporation have been moved from gog Shell Building, 
Houston, Texas, to 1041 Mellie Esperson Building. 


Louts StatuaM has left the Humble Oil and Refining Company to accept a position 
as Physicist at the University of California. His address is 921 Harbor View Drive, San 
Diego, California. 


H. C. Rercuert should be addressed c/o Isaac A. Manning, Barranquilla, Colom- 
bia, S. A. 


O. C. OLson has accepted a position with the Oil and Gas Division of the Texas 
Railroad Commission as a Petroleum Engineer. He is located at Brady, Texas. 


M. M. Stornick, former Editor of Gzopuysics under whom this department was 
started, has reported for duty at the Naval Ordnance Laboratory in Washington, D. C. 


S. C. StoneHam, Party Chief with Geophysical Service, Inc., is now located at 2407 
Jefferson Park Avenue, Charlottesville, Virginia. 


Maurice SKtar, of the Shell Oil Company, Inc., has been transferred from Big 
Lake to Box 566, Pearsall, Texas. 
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LAWRENCE A. GOEBEL, of the Carter Oil Company, has been transferred from De- 
catur, Illinois, to P. O. Box 548, Ferriday, Louisiana. 


WESLEY, J. RICKEL is acting as Observer with the Geophysical Research Corpora- 
tion, Box 2040, Tulsa, Oklahoma. 


Lr. LEonARD F. Unric is serving with the Maintenance Battalion, 5th A. D., at 
Camp Cooke, Lompoc, California. 


T. J. O'DONNELL of the Gulf Research and Development Company, is located at 
P. O. Box 733, Laurel, Mississippi. 


Cuar_es E. Houston has left the Humble Oil and Refining Company to accept 
a full-time position in the Department of Physics of the University of Houston, Hous- 
ton, Texas. 


S. ZIMERMAN, of the Carter Oil Company, has been transferred from Jackson, Mis- 
sissippi, to Box 441, Hope, Arkansas. 


THE SOCIETY OF EXPLORATION GEOPHYSICISTS 
CONSTITUTION AND BY-LAWS 
(As amended to April, 1942) 


I. NAME 


This Association shall be called the “Society of Exploration Geophysicists.” 


II. OByEctT 


The object of this Association is to promote the science of geophysics especially as 
it applies to exploration. 


III. MEMBERSHIP 


A—Members 


I. 


2. 


Any geophysicist of recognized standing shall be eligible to membership, whether 
or not he is engaged in petroleum geophysics or in geophysical prospecting. 

Any physicist, mathematician, geologist, or engineer of recognized standing 
who is investigating a geophysical problem or problems shall be eligible to 
membership. 


. Any geologist of recognized standing who is not a geophysicist of recognized 


standing, but who is acting as chief, or division chief in charge of geophysical 
surveying, shall be eligible to membership. 


B—Associates 

Any geophysicist, geologist, engineer, mathematician, or physicist who is a 
graduate of an institution of recognized scientific standing who is engaged in 
geophysics and who is not eligible to membership shall be eligible to associate 


membership. 


. Any person who is not a graduate of an institution of recognized scientific 


standing and who is not eligible to membership and who has shown ability in a 
routine phase of geophysical work and who is in a position of responsibility in 
geophysical surveying shall be eligible to associate membership. 


. Any representative or employee of a firm or corporation engaged in supplying 


equipment or supplies to the geophysical prospecting industry, provided such 
individual is the graduate of an institution of recognized scientific standing, or 
provided otherwise that such individual is qualified by adequate technical 
experience, shall be eligible to associate membership. 


. Associate members shall be known as associates. 
. Associates shall enjoy all the privileges of membership in the Society, save that 


they shall not hold office, sign applications for membership, or vote; neither 
shall they have the privilege of advertising their affiliation with the Society in 
professional cards or professional reports or otherwise. 


. The Executive Committee may advance to active membership, without the 


formality of application for such change, those associates who have, subsequent 
to election, fulfilled the requirements for active membership. 


C—Student Members 


I. 


Any graduate or undergraduate student in residence at a recognized university 
or college shall be eligible for student membership. 


2. Student membership shall terminate one year after a student member ceases to 


be enrolled at a recognized university or college. 
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3. The privileges of student members shall be the same as those of associate 
members, 


D—Election to Membership 

1. Every candidate for admission as a member or associate shall submit a formal 
application on an application form authorized by the Executive Committee, 
signed by him, and endorsed by not less than three members who are in good 
standing, stating his training and experience and such other facts as the Execu- 
tive Committee shall from time to time prescribe. Provided the Executive Com- 
mittee, after due consideration, shall judge that the applicant’s qualifications 
meet the requirements of the Constitution, they shall cause to be published in 
an appropriate publication the applicant’s name and the names of his sponsors. 
If, after at least thirty days have elapsed since such publication, no reason is 
presented why the applicant should not be admitted, he shall be deemed 
eligible to membership or to associate membership, as the case may be, and 
shall be notified of his election. Provided the Executive Committee may accept 
applicants for either class of membership, who, because of their geographical 
location or because of their lack of acquaintance with active members of this 
Society, are unable to secure qualified sponsors, but who otherwise are eligible 
for membership. Applicants for student membership may be accepted as student 
members upon recommendation of a faculty member of the school they are 
attending, and with the unanimous approval of the Executive Committee. 

. An applicant for membership, on being notified of his election in writing, shall 
pay full membership dues for the current year and on making such payment 
shall be entitled to receive the regular Society publications for that year. 
Unless payment of dues is made within thirty (30) days by those living within 
the continental United States and within ninety (go) days by those living 
elsewhere, after notice has been mailed of his election, the Executive Committee 
may rescind the election of the applicant. Upon payment of dues, each applicant 
for membership shall be furnished with a membership card for the current year 
and until such written notice and card are received, he shall in no way be 
considered a member of the Society. 

E—Honorary Members 
1. The Executive Committee may from time to time and by unanimous action elect 
as honorary members persons who have contributed distinguished service to 
the cause of geophysics. Honorary members shall not be required to pay dues. 


F—Code of Ethics 
1. Membership of any class shall be contingent upon conformance with the estab- 


lished principles of professional ethics. 


IV 


A—Resignation—Suspension—Expulsion 
1. Any member or associate may resign from the Society at any time. Such resig- 
nation shall be in writing and shall be accepted by the Executive Committee, 
subject to the payment of all outstanding dues and obligations of the resigning 
member or associate. 
. Any member or associate who is more than one year delinquent (in arrears) 
in payment of dues shall be suspended from the Society. Any delinquent or 
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suspended member or associate, at his own option, may request in writing that 
he be dropped from the Society and such request shall be granted by the Execu- 
tive Committee. Any member or associate more than two years in arrears shall 
be dropped from the Society. The time of payment of delinquent dues for either 
one year or two years may be extended by unanimous vote of the Executive 
Committee. 

3. Any member or associate who resigns or who is dropped under the provisions of 
Sections 1 and 2 of this article ceases to have any rights in the Society and 
ceases to incur further indebtedness to the Society. 

4. Any person who has ceased to be a member or associate under Section 1 or 
Section 2 of this article may be re-instated by unanimous vote of the Executive 
Committee subject to the payment of any outstanding dues and obligations 
which were incurred prior to the date when he ceased to be a member or associ- 
ate of the Society. 

5. Any member or associate who, after being granted a hearing by the Executive 
Committee, shall be found guilty of a violation of the code of ethics of this 
Society or shall be found guilty of a violation of the established principles of 
professional ethics or shall be found guilty of having made a false or misleading 
statement in his application for membership in the Society, shall be asked to 
resign from the Society by unanimous vote of the Executive Committee. The 
decision of the Executive Committee in all matters pertaining to the interpreta- 
tion and execution of the provisions of this section shall be final. 


V. OFFICERS AND THEIR DUTIES 
A—Officers 

1. The officers of the Society shall be a President, a Vice-President, a Secretary- 
Treasurer, and an Editor. These, together with the Past President, shall con- 
stitute the Executive Committee and managers of the Society. Officers shall be 
elected for a term of one year with the exception of the Editor, who shall be 
elected for a term of three years. 

2. A ticket of nominations for officers and representatives shall be prepared by 
a Committee on Nominations which shall consist of the President and the two 
qualified Past Presidents in order of precedence. Favorable action of this Com- 
mittee shall be two out of three votes. This Committee shall nominate two 
candidates for President, one of whom shall be the retiring Vice-President, and 
two candidates for each of the other offices to be filled. 

The ticket prepared by the Nominating Committee, as well as any other 
ticket presented in writing and signed by at least twenty members in good 
standing, shall be submitted to the President by November 30. The President 
shall secure written acceptance of nominations from each candidate. The names 
of nominees shall be placed on the official ballot which shall be mailed to each 
member not later than December 15. Names of candidates shall be listed 
alphabetically except that the name of the retiring Vice-President shall occupy 
first place cn the ticket. Blank spaces shall be provided on the ballot for writing 
in nominations for each office. 

Each member shall cast one vote for each officer to be elected and shall return 
his ballot to the Secretary-Treasurer. Ballots to be valid must be enclosed in an 
envelope carrying on the outside the written signature of the member submitting 
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the ballot, and must be received by the Secretary-Treasurer at his officially 
recognized address not later than ten days prior to the Annual Meeting. Only 
ballots received from members in good standing as of a date ten days prior to 
the Annual Meeting shall be valid. The Secretary-Treasurer shall indicate which 
ballots are valid by initialing the enclosing envelope and shall deliver the ballots 
unopened to the Nominating Committee for counting. 
The Nominating Committee shall count the ballots immediately before the 
Annual Business Meeting. The candidate receiving the greatest number of valid 
votes cast for an office shall be declared elected to that office. In case of a tie, 
the Nominating Committee shall decide by a secret vote which of the candidates 
tied shall be elected. 
The announcement of the election shall be made as the first of the new busi- 
ness considered at the Annual Meeting. 

3. No one shall hold the office of President for two consecutive years and no one 
shall hold any other office for more than two consecutive years except that of 
Editor. 


B—Duties of Officers 

1. The President shall be the presiding officer at all the meetings of the Society, 
shall take cognizance of the acts of the Society and of its officers, shall appoint 
such committees as are required for the purposes of the Society, and shall dele- 
gate members to represent the Society. He may, at his option, serve on, and 
may be chairman of, any committee. 

2. The Vice-President shall assume the office of President in case of a vacancy 
from any cause in that office and shall assume the duties of President in case 
of the absence or disability of the latter. He shall also act as Chairman of the 
Program and Arrangements Committee of the Society. 

3. The Secretary-Treasurer shall assume the duties of the President in case of the 
absence of both the President and Vice-President. He shall have charge of the 
financial affairs of the Society and shall annually submit reports as Secretary- 
Treasurer covering the fiscal year, which he shall arrange to have published in 
the next regular issue of the recognized organ of the Society. He shall receive all 
funds of the Society and, under the direction of the Executive Committee, shall 
disburse all funds of the Society. He shall cause an audit to be prepared an- 
nually by a public accountant at the expense of the Society. He shall give a 
bond, and shall cause to be bonded, all employees to whom authority may be 
delegated to handle Society funds. The amount of such bonds shall be set by the 
Executive Committee and the expense shall be borne by the Society. The funds 
of the Society shall be disbursed by check as authorized by the Executive Com- 
mittee. 

4. The Editor shall be in charge of the editorial business, shall submit an annual 
report of such business, shall have authority to solicit papers and material for 
the regular Society publication and for special publications, and may accept 
or reject material offered for publication. He may appoint associate, regional, 
and special editors. 

5. In the event of absence, disability or resignation of any member of the Executive 
Committee, the vacant place in the Committee shall be automatically filled by 
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the next preceding qualified Past President who shall be known as Prior Past 
President. 

6. The officers shall assume the duties of their respective offices immediately after 
the Annual Meeting following their election. 


VI. EXEcUTIVE COMMITTEE—MEETINGS AND DUTIES 


A—Executive Committee 
1. The Executive Committee shall consist of the President, Past President, Vice- 


President, Secretary-Treasurer, and Editor. 


B—Meetings and Duties 

1. The Executive Committee shall meet preceding the annual meeting and at the 
call of the President may hold meetings when and where thought advisable, 
to conduct the affairs of the Society. A joint meeting of the outgoing and 
incoming Executive Committees shall be held immediately after the close of the 
annual Society business meeting. Members of the Executive Committee may 
vote by proxy, mail or person on all matters which cal] for favorable or unani- 
mous action. 

. The Executive Committee shall consider all nominations for membership and 
pass on the qualifications of the applicants, shall have control and management 
of the affairs and funds of the Society, shall determine the manner of publica- 
tion, and shall designate the place of the annual meeting. They are empowered 
to establish a business headquarters for the Society, and to employ such persons 
as are needed to conduct the business of the Society. They are empowered to 


make investments of both general and special funds of the Society. Trust funds 
may be created giving to the trustees appointed for such purpose such discretion 
as to investments as seems desirable to the Executive Committee to accom- 
plish any of its objects and purposes, but no such trust funds shall be created 
unless they are revocable upon ninety (go) days’ notice. Favorable action shall 
consist of a favoring three out of five votes. Committee members may vote by 
person, proxy or in writing. 


Vil 
A—Meetings 

1. The Society shall hold at least one stated meeting each year, which shall be the 
annual meeting. This meeting shall be held at a time and place designated by 
the Executive Committee. At this meeting the proceedings of the preceding 
meeting shall be read, Society business shall be transacted, scientific papers 
shall be read and discussed, and results announced of the mail ballot for officers 
for the ensuing year. 


VIII 
A—Amendments 
1. Amendments to this Constitution may be proposed by a resolution of the Execu- 
tive Committee, by a Constitutional Committee appointed by the President, 
or in writing by any ten members of the Society. All such resolutions or proposals 
must be submitted at the annual meeting of the Business Committee of the So- 
ciety as provided in the By-Laws, and only the Business Committee shall make 
recommendations concerning proposed constitutional changes at the annual 
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Society business meeting. If such recommendations by the Business Committee 
shall be favorably acted on at the annual Society business meeting, the Secre- 
tary-Treasurer shall arrange for a ballot of the membership by mail within thirty 
(30) days after said annual Society business meeting, and a majority vote of the 
ballots received within ninety (90) days of their mailing shall be sufficient to 
amend. The legality of all amendments must be determined by the Executive 
Committee prior to balloting. 


B—Parliamentary Regulations 
1. The rules contained in ‘‘Robert’s Rules of Order, Revised” shall govern the 
Society in all cases to which they are applicable, and in which they are not in- 
consistent with the By-Laws or the special rules of order of this Society. 


BY-LAWS 


I 
A—Dues 

1. The fiscal year of the Society shall correspond with the calendar year. 

2. The annual dues of members of the Society shall be five dollars ($5.00). The 
annual dues of associates for not to exceed six years after election shall be four 
dollars ($4.00); thereafter, the annual dues of such associates shall be five dollars 
($5.00). The annual dues of student members shall be three dollars ($3.00). The 
annual dues are payable in advance on the first day of each calendar year. A 
bill shall be mailed to each member, associate, and student member before 
January first of each year, stating the amount of the annual dues and the 


penalty and conditions for default in payment. Members, associates, or student 
members who shall fail to pay their annual dues by March fifteenth shall not 
receive further copies of the regular Society publication, shall not be eligible to 
vote, nor shall they be privileged to buy Society special publications at prices 
made to the membership, until such arrears are met. 


II 
A—Publications 
1. The proceedings of the annual meeting, the Constitution and By-Laws, and the 
papers presented at such meeting shall be published at the discretion of the 
Executive Committee in the regular publication or in such other form as the 
Executive Committee may decide best meets the needs of the membership of 
the Society. 

. The payment of annual dues for any fiscal year entitles the member or associate 
to receive without further charge the regular publication of the Society for that 
year. $3.00 ($3.50 foreign) of the annual dues is set aside for pepnene of the an- 
nual subscription to the Society Journal, GEopxHysics. . 

. The Executive Committee may authorize the printing of special publications 
to be financed by the Society from its general, publication or special funds and 
offered for sale to members and associates in good standing at not less than the 
cost of publication and distribution. 


Ill 


A—Regional Sections, Technical Divisions, and Affiliated Societies 
1. Regional sections of the Society may be established provided the members of 
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such sections are members of the Society and shall perfect an organization and 
make application to the Executive Committee. The Executive Committee shall 
submit the application for action first by the Business Committee and second 
to the Society for a letter ballot as for a constitutional amendment, and pro- 
vided that the Society may revoke the charter of any regional section by similar 
action. 

2. The Executive Committee may arrange for the affiliation with the Society of 
duly organized groups or societies, which by object, aims, constitution, by-laws, 
or practice are developing the study of geophysics or petroleum technology. 
In like manner, the Executive Committee may arrange conditions for dissolu- 
tion of such affiliations. Affiliation with the Society need not prevent affiliation 
with other scientific societies. Members of affiliated societies who are not mem- 
bers of the Society, shall not have the privilege of advertising their affiliation 
with the Society on professional cards or otherwise, but shall be entitled to re- 
ceive the regular publication of the Society upon payment of a fee equivalent to 
that charged an associate member. 


IV 

A—District Representatives 
1. The Executive Committee may cause to be elected district representatives from 
districts which it shall define by a local geographic grouping of the membership. 
Such districts shall be redesignated and redefined by the Executive Committee 
as often as seems advisable. Each district shall be entitled to one representative 
for each seventy-five members, but this shall not deprive any designated district 
of at least one representative. The representatives so apportioned shall be chosen 
from the membership of the district by a written ballot arranged by the Execu- 
tive Committee. They shall hold office for two years, their term of office expiring 
at the close of the annual meeting. In so far as possible, these two-year terms 

shall be overlapping. 


V 
A—Business Committee 

1. The Executive Committee may institute a Business Committee to act as a coun- 
cil and advisory board to the Executive Committee and the Society. This 
Committee shall consist of the Executive Committee, the Prior Past President, 
and the district representatives. Favorable action by the Business Committee 
shall consist of a favorable two-thirds vote. The Secretary-Treasurer shall act as 
Secretary of the business meeting. Members of this Committee may vote in per- 
son, by proxy, or in writing. The Business Committee shall meet the day before 
the annual meeting at which time proposed changes in the Constitution or By- 
Laws may be considered, all old and new business shall be discussed, and recom- 
mendations may be voted for presentation at the annual meeting. 


VI 
A—Amendments 
1. These by-laws may be amended by favorable vote of the Business Committee, 
providing that such changes have been recommended and legality approved by 
the Executive Committee. Such amendments shall be published in the next 
regular number of the recognized organ of the Society. 
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